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Abstract 

Monostatic  radars  employ  pulsed  waveforms  to  achieve  receiver /transmitter 
isolation.  Pulsed  waveforms,  however,  are  intrinsically  ambiguous  in  range  and 
Doppler.  The  mitigation  of  these  inherent  ambiguities  has  been  the  focus  of  radar 
waveform  designers  for  many  decades.  In  1962  Palermo  used  two  conjugate  linear 
frequency  modulated  (LEM)  pulses  to  demonstrate  a  concept  that  reduces  ambigu¬ 
ous  energy  in  processed  radar  return  signals.  This  approach,  known  as  Non-linear 
Ambiguity  Suppression  (NLS),  acts  on  the  radar  return  for  a  pulse-coded  waveform 
by  compressing  undesirable  pulses  and  suppressing  them  using  saturation  threshold¬ 
ing.  This  process  effectively  reduces  the  energy  contribution  of  undesirable  pulses 
in  the  resultant  radar  return.  In  this  way  Palermo  was  able  to  provide  a  2-fold  im¬ 
provement  in  range  ambiguities,  however  the  use  of  conjugate  LEM  pulses  does  not 
extend  to  larger  symbol  families,  thus  having  severe  limitations. 

In  order  to  overcome  the  limitation  on  symbol  family  size,  NLS  radar  waveforms 
can  comprise  pseudo-random  discrete  codes,  such  as  Gold  codes,  that  exhibit  strong 
mutual  dispersion  properties.  However,  no  theorem  exists  that  directly  tackles  the 
problem  of  designing  optimal  mutually  dispersive  symbols.  This  thesis  analyzes  the 
Brown  theorem,  a  recent  theorem  devised  in  support  of  NLS.  The  Brown  theorem 
utilizes  a  rms  time  duration  metric  to  design,  in  the  frequency  domain,  optimal 
mutually  dispersive  symbols  for  any  sized  pulse  family.  The  Fourier  transform  of 
the  rms  time  duration  is  considered  because  it  neatly  isolates  the  envelope  and 
phase  components  of  the  NLS  waveform  thereby  permitting  the  systematic  design  of 
mutually  dispersive  pulses.  Optimality  in  compression  is  found  from  the  solution  to 
a  constrained  optimization  problem,  whilst  mutuality  in  dispersion  is  ensured  via  a 
prudent  selection  of  phase-rate  functions. 


xiii 


This  thesis  signihcantly  advances  research  towards  the  implementation  of  op¬ 
timal  NLS  waveforms  by  analyzing  the  Brown  theorem.  The  Brown  theorem  is 
reintroduced  with  the  use  of  simplihed  linear  algebraic  notation.  A  methodology 
for  Brown  symbol  design  and  digitization  is  provided,  and  the  concept  of  disper¬ 
sive  gain  is  introduced.  Numerical  methods  are  utilized  to  design,  synthesize,  and 
analyze  Brown  symbol  performance.  The  theoretical  performance  in  compression 
and  dispersion  of  Brown  symbols  is  demonstrated  and  is  shown  to  exhibit  signihcant 
improvement  compared  to  discrete  codes.  As  a  result  of  this  research  a  process  is 
presented  for  the  design  of  optimal  mutually  dispersive  symbols  for  any  sized  family. 
In  other  words,  the  limitations  imposed  by  conjugate  LFM  are  overcome  using  NLS 
waveforms  that  provide  an  effective  A^-fold  increase  in  radar  unambiguous  range. 
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AN  ANALYSIS  OF  MUTUALLY  DISPERSIVE  BROWN 
SYMBOLS  FOR  NON-LINEAR  AMBIGUITY  SUPPRESSION 


/.  Introduction 

1.1  Research  Motivation 

The  radar  system  is  primarily  employed  for  providing  information  about  the 
presence  of  a  target,  and  its  associated  parameters  such  as  size,  location  and  motion. 
This  information  is  sought  after  in  a  hierarchy  from  target  detection  to  parameter 
estimation  with  waveform  resolution  being  the  governing  metric  of  precision  [31]. 
Accordingly,  radar  waveforms  are  chosen  based  on  application  so  that  they  deliver 
sufficient  precision  in  resolution.  With  their  origins  in  continuous  wave  systems, 
radars  quickly  evolved  into  pulsed  systems  before  the  end  of  World  War  II.  Monos¬ 
tatic  radars  (collocated  transmitter  and  receiver)  utilize  pulsed  waveforms  so  that  re¬ 
turns  are  received  free  of  interference  from  radar  transmissions.  Originally  extremely 
narrow  pulses  were  used  to  achieve  high  resolution,  however,  due  to  hardware  satu¬ 
ration  limitations  this  resulted  in  low  transmit  power  and  was  deemed  ineffective.  In 
1960,  Klauder  et  ah  [17]  introduced  the  linear  frequency  modulated  (LFM),  or  chirp, 
pulse  which  exploited  spectral  content  to  achieve  pulse  compression;  this  marked  the 
beginning  of  the  era  of  pulse  compression  waveform  design. 

A  pulsed  radar  transmits  short  pulses  of  electromagnetic  energy  at  a  rate  known 
as  the  Pulse  Repetition  Frequency  (PRF).  This  rate  varies  based  on  operational  re¬ 
quirements  but  intrinsically  leads  to  ambiguities  in  range  (time  delay)  and  Doppler 
(velocity)  estimation.  By  changing  the  PRF  these  ambiguities  may  be  controlled  but 
not  necessarily  removed.  A  number  of  different  methods  exist  for  reducing  the  de¬ 
gree  of  ambiguity  in  radar  returns.  However,  the  previously  transmitted  pulse  that 
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generated  the  return  under  consideration  must  be  determined  to  entirely  remove 
range  ambiguities.  The  ability  to  identify  a  particular  member  of  a  set  of  uniquely 
coded  pulses  from  all  others  in  the  set  is  known  as  pulse  discrimination.  This  proves 
invaluable  in  many  applications,  such  as  pulsed-Doppler  radar  and  wireless  spread- 
spectrum  communications  [30].  Pulse  compression  and  pulse  discrimination  are  re¬ 
lated  through  mutual  dispersion.  A  set  of  symbols  having  crosscorrelation  functions 
with  equivalent  dispersive  characteristics  are  said  to  be  mutually  dispersive.  If  the 
same  set  of  symbols  exhibit  good  pulse  compression,  i.e.,  have  focused  autocorrela¬ 
tion  functions,  then  pulse  discrimination  becomes  rudimentary.  In  the  ideal  case  the 
matched  output,  or  autocorrelation  function,  would  be  the  dirac  delta  function  (all 
energy  concentrated  at  a  single  point).  Clearly  this  is  not  realizable,  therefore  highly 
dispersive  crosscorrelation  functions  are  required  to  enable  pulse  discrimination  using 
a  hlter  and  threshold  detector. 

Intent  on  improving  unambiguous  range,  Palermo  [28]  discovered  a  means  of 
combining  pulse  compression  and  pulse  discrimination  into  one  technique,  known  as 
Non-linear  Suppression  (NLS).  NLS  utilizes  pulse  discrimination  to  reduce  ambigu¬ 
ous  energy  in  a  radar  return  with  the  added  beneht  of  exploiting  pulse  compression 
to  achieve  high  range  resolution.  Each  channel  of  a  NLS  system  is  assigned  a  unique 
symbol;  the  radar  return  signal  is  passed  through  the  channel  producing  an  output 
indicating  the  presence  of  that  symbol  alone.  Consequently,  a  series  of  suppression 
operators  are  employed  in  each  channel  to  sequentially  suppress  the  response  of  all 
undesirable  symbols.  These  suppression  operators  are  the  kernel  of  NLS  that  exploit 
the  pulse  discrimination  concept  to  suppress  undesirable  symbol  responses. 

Given  a  suitable  mutually  dispersive  symbol  set,  Palermo’s  suppression  kernel 
comprises  three  components:  pre-hlter,  non-linearity,  and  post-hlter.  The  pre-hlter 
is  matched  for  a  particular  undesirable  symbol,  resulting  in  a  highly  focused  autocor¬ 
relation  response,  and  dispersed  crosscorrelation  response  for  all  other  symbols.  The 
pre-hltered  signal  then  undergoes  a  saturating  non-linearity  that  clips  (limits)  the 
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input  above  a  prescribed  threshold.  This  effectively  removes  a  good  deal  of  undesir¬ 
able  symbol  energy  without  affecting  other  symbols.  Finally,  the  suppressed  signal  is 
passed  through  a  post-filter  that  inverts  the  pre-filter  effect.  Sequentially  performing 
this  NFS  operation  for  all  but  the  desired  pulse  reduces  the  ambiguous  energy  in  the 
resultant  signal.  When  the  desired  symbol  is  finally  focused,  or  matched-filtered, 
the  ambiguous  energy  present  in  the  processed  radar  return  is  significantly  reduced. 
Therefore,  NFS  is  an  interpulse  coding  scheme  that  exploits  pulse  discrimination  to 
provide  a  means  for  reducing  pulsed  radar  range  ambiguities. 

A  fundamental  prerequisite  to  employing  NFS  is  having  a  set  of  mutually 
dispersive  symbols.  In  [5]  Brown^  presents  a  theorem  on  optimal  mutually  dispersive 
symbols.  As  stated  earlier,  mutually  dispersive  symbols  may  be  characterized  by  the 
dispersive  characteristics  of  their  correlation  functions.  By  considering  the  Fourier 
transform  of  the  rms  time  duration  of  the  correlation  functions.  Brown  devises  a 
theory  for  the  spectral  design  of  optimal  mutually  dispersive  symbols.  In  this  fashion, 
the  rms  time  duration  neatly  isolates  the  envelope  and  phase  components  of  the 
correlation  function  thereby  permitting  the  systematic  design  of  mutually  dispersive 
pulses.  This  formulation  indicates  that  the  optimal  spectral  envelope  is  the  solution 
to  a  constrained  optimization  problem,  and  the  optimal  phase-rate  functions  are 
linear  combinations  of  a  weighted-orthonormal  basis.  The  search  for  the  optimal 
linear  weighting  of  the  basis  functions  spawns  a  problem,  coined  by  Brown  as  the 
Hermit  Problem.  Brown  states  that  the  crosscorrelation  dispersion  of  such  symbols 
is  independent  of  basis  and  proportional  to  the  number  of  symbols,  N,  in  the  set  via 
2N/{N  —  1).  However,  the  theory  for  the  design  of  an  optimal  mutually  dispersive 
set  of  symbols,  known  as  Brown  symbols,  remains  untested. 

1.2  Research  Goal 

The  goal  of  this  research  is  to: 
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1.  Explain  the  formulation  of  Brown’s  theorem  [5]  whilst  introducing  a  simplihed 
linear  algebraic  notation; 

2.  Describe  a  systematic  process  for  designing  Brown  symbols  through  an  analysis 
of  Brown’s  work  and  by  theoretical  extension; 

3.  Provide  a  theoretical  analysis  for  the  digitization  of  Brown  symbols; 

4.  Evaluate  performance  of  Brown  symbols,  comprising  dispersive  gain,  through 
simulation  and  numerical  approximation;  and 

5.  Compare  the  performance  of  various  Brown  symbols  with  discrete  codes  through 
simulation. 

1.3  Thesis  Organization 

Chapter  II  introduces  background  theory  and  presents  a  review  of  current  lit¬ 
erature.  Due  to  the  theoretical  nature  of  this  thesis,  a  functional  analysis  overview  is 
provided  dehning  nomenclature  and  metrics.  Numerical  approximation  techniques 
and  the  role  of  the  Discrete  Eourier  Transform  (DFT)  are  discussed.  The  reference 
radar  model  and  generic  form  for  the  radar  return  signal  are  dehned.  Bandpass  sys¬ 
tems  and  analytic  notation  are  discussed,  including  consideration  of  optimal  hltering 
in  white  noise.  Pulsed  waveform  design  is  discussed  in  light  of  range  and  Doppler  am¬ 
biguities  of  monostatic  pulse  Doppler  radar,  and  metrics  for  signal  characterization 
are  dehned,  including  rms  time  duration.  Finally,  ambiguity  suppression  is  reviewed, 
consideration  is  given  to  PRF  techniques,  and  the  NLS  concept  is  explained. 

Chapter  III  explains  Brown’s  theorem  and  introduces  the  symbol  design  method¬ 
ology.  A  brief  overview  of  the  NLS  concept  is  provided.  The  scope  of  the  research  is 
outlined  as  relating  to  the  approach  and  associated  synthesis  constraints.  The  radar 
model  is  related  to  NLS  by  analytically  dehning,  as  in  [5] ,  the  problem  of  designing  a 
symbol  set  suited  to  NLS,  and  the  Brown  NLS  system  is  illustrated.  Brown’s  theorem 
is  derived  by  tying  together  the  optimal  envelope,  phase-rate  functions.  Hermit  Prob- 
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lem  solutions,  dispersive  gain  and  the  theoretically  optimal  performance.  Notation 
for  biphase  codes  is  dehned  relating  the  suitability  of  these  codes  to  NLS.  Finally, 
metrics  used  to  characterize  and  contrast  Brown  symbols  are  discussed,  providing 
an  overview  of  the  numerical  approximation  approach. 

In  Chapter  IV,  Brown  symbols  are  devised  analytically  and  then  constructed, 
providing  results  from  Matlab®  simulations.  The  design  requirements  are  consid¬ 
ered  for  each  Brown  symbol  components:  the  Hermit  Problem  solution,  phase-rate 
basis  function  design,  phase  function  calculation,  and  symbol  synthesis.  Results  are 
provided  demonstrating  the  sensitivity  of  Brown  symbols  to  the  choice  of  phase-rate 
basis  functions  and  envelope;  physical  realizability  is  also  discussed.  Finally,  the 
performance  of  Brown  symbols  is  contrasted  with  that  of  biphase  symbols  suited  to 
NLS. 

Chapter  V  provides  a  summary  of  the  signihcant  research  hndings.  The  re¬ 
search  goal  is  restated  and  the  conclusive  results  summarized.  The  signihcance  of 
this  research  is  discussed  and  recommendations  for  future  research  provided. 
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II.  Background 


2. 1  Introduction 

This  chapter  introduces  background  theory  and  presents  a  review  of  current 
literature.  Due  to  the  theoretical  nature  of  this  thesis  Section  2.2  provides  a  func¬ 
tional  analysis  overview,  dehnes  nomenclature  and  metrics,  and  discusses  Numerical 
approximation  techniques  and  the  role  of  the  Discrete  Fourier  Transform.  Section  2.3 
dehnes  the  reference  radar  model  and  generic  form  for  the  radar  return  signal,  dis¬ 
cusses  bandpass  systems  and  analytic  notation  and  considers  optimal  hltering  in 
white  noise.  Section  2.3.3  discusses  pulsed  waveform  design  in  light  of  range  and 
Doppler  ambiguities  of  monostatic  pulse  Doppler  radar,  and  dehnes  signal  associated 
metrics  for  characterizing  radar  signals.  Finally,  Section  2.4  reviews  ambiguity  sup¬ 
pression,  considering  Pulse  Repetition  Frequency  (PRF)  techniques  and  Non-linear 
Suppression  (NFS). 

2.2  Signal  and  Waveform  Analysis 

2.2.1  Signal  Spaces.  The  construction  and  characterization  of  signal  spaces 
is  fundamental  to  this  thesis.  A  signal  is  dehned  as  “a  quantity  which,  in  some 
manner,  conveys  information  about  the  state  of  a  physical  system  [11:1].”  A  signal 
space,  therefore,  is  the  space  in  which  a  given  signal  can  be  considered  a  single 
point  or  entity.  The  signals  considered  in  this  thesis  are  associated  with  the  radar 
system  depicted  in  Figure  2.2,  and  are  denoted  as  functions  of  either  time,  x(t), 
or  frequency,  X{f).  With  the  constraint  of  physical  realizability  enforced,  these 
signals  are  constrained  to  exist  within  the  union  of  the  sets  of  energy-limited  and 
band-limited  signals,  that  is  {x  :  a;  G  Se{K)  [11:5]  where. 

Energy-limited  signals; 
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Band-limited  signals; 


/oo  ^ 

x(t)e-’^'>'dt  =  0  for  all  |/|  >  - 


(2,2) 


The  populous  of  physically  realizable  signals  is  considerable  and  therefore  re¬ 
quires  a  means  of  characterizing  the  different  elements  for  comparison.  Any  func¬ 
tional  that  maps  elements  or  pairs  of  elements  onto  the  real  line  provides  a  means 
of  characterization.  A  measure,  dehned  as  distance  [11:20],  is  used  to  produce  a  real 
positive  number.  Provided  that  the  functional,  d  :  {x,  y}  TZ,  exhibits  the  normal 
properties  of  distance  then  it  is  commonly  known  as  a  metric,  and  when  combined 
with  a  set  of  elements  results  in  a  metric  space  [11]. 

Of  all  possible  metrics,  the  analysis  undertaken  in  this  thesis  requires  only  an 
understanding  of  the  Euclidean  metric,  norm,  and  inner  product.  When  considering 
vectors  in  A^-dimensional  space,  or  TZ^ ,  with  vectors  x  =  [xq  xi  ...  xtv-i]^  and 
y  =  [yo  yi  •  •  •  yw-i]"^,  the  Euclidean  metric  (i2(x,y)  is  defined  as  [11] 


c^2(x,y)  =  ||x-y|| 

■  N-l  1  5 

.  k=0 


(2.3) 


When  considering  the  set  of  real  or  complex  time  functions,  x{t)  and  y{t), 
dehned  over  the  interval  {t  \  a  <t  <h},  the  Euclidean  metric  becomes  [11] 


d2{x,y) 


1 

2 


x(t)  —  y(t)\‘^dt 


(2.4) 
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The  norm  determines  the  “size”  of  a  particular  element  in  a  set  and  therefore 
is  a  particularly  useful  variant  of  the  Euclidean  metric.  This  metric  (i2(x,0)  =  ||x|| 
is  dehned  for  real  or  complex  vectors  as 


N-l  1  2 

A:=0 

1 

and  for  arbitrary  time  functions  as 
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x(t)\‘^dt 


(2.5) 


(2.6) 


Any  complex  or  real  time  function  with  hnite  size  belongs  to  the  set  of  hnite 
energy  signals,  Se{K),  dehned  in  Eq  (2.1).  To  complete  the  discussion  of  signal 
spaces,  a  metric  describing  the  relationship  between  pairs  of  vectors  is  dehned  in  the 
form  of  an  inner  product.  The  inner  product  of  a  pair  of  real  or  complex  vectors, 
and  arbitrary  time  functions  is  dehned  in  Eq  (2.7)  and  Eq  (2.8),  respectively: 


(x,y) 

{x,y) 


N-l 

k=0 


x(t)  y*(t)dt. 


(2.7) 

(2.8) 


In  some  applications  it  may  be  benehcial  to  emphasize  a  certain  time,  or  frequency, 
window  in  calculating  element  pair  relationships.  The  weighted  inner  product  is 
designed  to  provide  this  utility  and  is  described  by  [11] 


{x,y)w=  f  w{t)x{t)y*{t)dt.  (2.9) 

J  a 


The  norm  metric,  described  earlier,  in  conjunction  with  a  linear  space  forms 
a  normed  linear  space  with  ||x||  =  (x,x)5.  Also,  describing  an  inner  product  in  this 
way  leads  to  the  concept  of  orthogonality.  The  inner  product  between  two  vectors  is 
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related  to  the  norm  squared  of  the  vector  x+y  as  follows: 


y||^  =  ||x||^  +  ||yr  +  2xV 


(2.10) 


The  vector  x+y  is  a  vector  completing  a  triangle  having  sides  x  and  y.  The 
Pythagorean  Theorem  states  that  if  the  square  of  the  hypotenuse  equals  the  sum  of 
the  square  of  the  sides,  then  the  sides  are  orthogonal.  From  the  equality  described 
above,  this  implies  that  if  x^y  =  0,  then  x  is  orthogonal  to  y. 


Not  surprisingly,  a  signal  space  may  be  represented  by  far  fewer  elements  than 
are  contained  in  the  entire  space.  The  signal  space  described  by  all  linear  combina¬ 
tions  of  any  N  elements  is  called  a  linear  subspace,  denoted  S'Ar(x)  and  dehned 
as  [11] 


*S'v(x) 


N-l 

fc=0 


(2.11) 


Any  set  of  elements  from  a  signal  space  are  said  to  be  linearly  independent  if 
the  equality  =  0  satished  when  c^  =  0  for  all  k,  i.e.,  any 

cannot  be  written  as  a  linear  combination  of  the  other  elements.  When  this  is  the 
case,  {ufc}  is  known  as  a  basis  for  the  iV-dimensional  subspace  S'7v(x),  and  S'Ar(x) 
is  said  to  be  spanned  by  {u^}.  Therefore,  any  iV-linearly  independent  elements  in 
S'7v(x)  can  form  a  basis,  implying  that  a  basis  is  by  no  means  unique. 

Any  vector  by  this  construction,  x  G  5'Ar(x),  can  be  represented  by  a  linear 
combination  of  basis  elements.  The  weighting  of  the  l^^  basis  element  is  determined 
by  taking  the  inner  product  of  x  with  a  vector  that  is  orthogonal  to  all  other  basis 
elements,  and  that  normalizes  the  l^^  basis  element  as  detailed  below. 


N-l 

Cl  =  (x,Vi)  =  y^Cfc(Ufc,  V;) 
k=0 


(2.12) 
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This  process  is  simplified  greatly  by  choosing  an  orthonormal  basis  {q*;}.  That 
is  :  (qfc,q;)  =  Ski},  where  Ski  is  the  Kronecker  delta,  which  equals  unity  only 
when  k  =  I  and  zero  otherwise.  In  other  words,  q^  has  unit  size  and  is  orthogonal 
to  all  other  elements  in  the  basis.  In  this  way,  the  weighting  of  the  basis  element 
is  determined  by  simply  taking  the  inner  product  of  x  with  the  basis  element. 

A  basis  may  be  orthonormalized  using  the  Gram-Schmidt  orthonormalization 
procedure.  Given  a  set  of  iV-linearly  independent  vectors  {u^},  the  orthonormal  set 
{qj  is  calculated  by  normalizing  the  set  {w^}  obtained  by  [18:279]: 


Wo 

=  Uo 

Wi 

=  Ui-(ui,qo)qo 

W2 

=  U2  -  (U2,qi)qi - 

fc-i 

Wfc 

=  Ufc  -  ^(Ufc,q;)q; 

z=o 

where 

_  Wfc 

llWfcll 

(2.13) 


(2.14) 


2.2.2  Integral  Transforms.  Section  2.2.1  discussed  signal  spaces  and  met¬ 
rics  that  are  commonly  used  to  characterize  signals  that  exist  as  elements  within 
them.  This  section  explores  the  relationship  of  these  concepts  to  Integral  trans¬ 
forms. 

As  discussed  in  section  2.2.1,  linear  subspaces  may  be  described  by  a  linear 
combination  of  a  set  of  basis  elements  as  illustrated  in  Eq  (2.11).  The  study  of 
continuous  time  functions,  in  Se{K),  can  be  considered  analogous  by  viewing  an  N- 
dimensional  vector  as  a  discrete  approximation  to  a  continuous  time  function  whose 
approximation  error  decreases  as  the  sampling  interval  approaches  zero.  In  this  way 
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Eq  (2.11)  can  be  redefined  for  continnous  signals  by  the  fnnctional 


x{t)  =  (2-15) 

where  S  is  the  range  of  Note  that  c{^)  is  analogous  to  Ck  in  Eq  (2.11)  and  that 
called  the  basis  kernel,  is  the  time  continuous  basis  function  parameterized 
by  In  order  to  complete  this  analogy  c(.^)  is  described  as 


c(0  =  y 


(2.16) 


where  T  is  the  range  of  t.  This  equation  is  the  inner  product  of  x{t)  with 
known  as  the  reciprocal  basis  kernel.  For  linear  subspaces  this  is  the  set  of  basis 
elements  that  enables  the  recovery  of  the  original  basis  weighting  due  to  the  imposed 
orthogonality  relationship. 

When  a  complex  sinusoid  is  used  as  the  basis  kernel  the  well  known  Fourier 
transform  with  the  Fourier  transform  pair  is  described  here: 


T{xit)}  :  Xif)  ^ 
X-^{X{f)}  :  x{t)  t 


x{t)e 


(2.17) 

(2.18) 


The  frequency  content  in  a  continuous  time  function  is  divined  from  this  transform. 

There  are,  however,  properties  of  the  Fourier  transform  that  should  be  dis¬ 
cussed  in  order  to  relate  this  concept  to  an  analysis  of  a  radar  system.  For  x{t)  a 
real- valued  function,  the  complex  conjugate  of  Eq  (2.17)  shows  that  X*{f)  =  X(— /). 
This  is  significant  as  it  illustrates  a  symmetry  in  the  frequency  content  of  real  time 
functions,  and  therefore  a  redundancy.  That  is,  only  one-half  of  the  spectrum  is 
required  to  define  a  real-valued  signal.  Consequently,  an  analysis  of  real  signals  is 
simplified  by  representing  them  with  a  one-sided  spectral  signal,  or  analytic  signal. 
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One  approach  is  to  define  a  complex  valued  signal,  'ip(t)  G  C,  that  has  a  one¬ 
sided  spectrum  equivalent  to  that  of  the  real- valued  signal,  tI/(/)  =  X{f)  for  /  >  0, 
and  whose  real  part  is  equal  to  the  signal,  =  x(t).  For  example: 

ip{t)  =  x{t)  +  jx{t)  (2.19) 

^{f)=X{f)+jX{f).  (2.20) 


'!/(/)  is  to  be  one-sided,  and  has  arbitrarily  been  assumed  non-negative  for  positive 
frequencies.  Therefore,  the  spectrum  of  the  imaginary  part  of  the  analytic  signal, 
X(/),  must  cancel  the  spectrum  of  X{f)  for  negative  frequencies  and  as  such  is 
defined  as 


^(/) 


-jX{f)  />0 
jX{f)  /<0 


-jX{f)  Sgn(/), 


(2.21) 


where  sgn(/)  is  1  for  /  >  0  and  —1  for  /  <  0.  Therefore,  from  Eq  (2.20)  and 
Eq  (2.21)  the  spectral  density  of  the  analytic  signal  becomes  one-sided  and  is  of  the 
form 


^(/) 


2X(/)  />0 
0  /<0 


(2.22) 


The  time  function,  x{t),  can  be  determined  by  taking  the  inverse  Fourier  transform  of 
X{f)  defined  in  Eq  (2.21).  For  jsgn(/)}  =  l/7rf,  by  the  convolution  property 

of  the  Fourier  transform,  x{t)  is  given  by 


H{x(t)}  :  x(t)  =  x{t)'k  — 

_  I  r  x{s) 

J—OC  f  ^ 


(2.23) 


Eq  (2.23)  introduces  the  Hilbert  transform  definition,  another  integral  transform  as 
defined  in  Eq  (2.15).  In  this  case,  however,  l/(t  —  s)  is  the  basis  kernel  and  exhibits 
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a  singularity  at  s  =  t,  therefore,  define 


lim 

£^0+ 


2.2.3  Discrete  Signal  Approximations.  Although  much  of  the  later  dis¬ 
cussions  in  this  thesis  are  related  to  an  optimization  approach  regarding  continuous 
time  signals,  a  physical  implementation  is  likely  to  require  some  manner  of  digi¬ 
tization.  Therefore,  this  section  provides  the  digital  signal  processing  background 
necessary  to  facilitate  this  transition  and  to  introduce  some  numerical  approximation 
fundamentals. 

The  use  of  a  complex  exponential  as  a  basis  kernel,  as  seen  in  the  Fourier 
transform  introduced  in  Section  2.2.2,  is  commonly  considered  for  various  forms  of 
discrete  signals.  Signals  are  typically  taken  to  be  either  continuous  or  sampled  ap¬ 
proximations,  therefore  an  appropriate  revision  of  the  Fourier  transform  is  required. 
The  transform  pair  for  continuous  signals  in  time  and  frequency  is  described  in 
Eq  (2.17);  signals  that  are  continuous  and  periodic  in  time  may  be  represented  by  a 
Fourier  series  [25],  i.e.,  having  discrete  or  sampled  frequency  spectra.  The  duality  of 
frequency  and  time  signals,  as  related  via  the  Fourier  transform,  ensures  the  spectra 
of  a  time  sampled  signal  is  continuous  and  periodic;  being  related  via  the  Discrete- 
Time  Fourier  Transform  (DTFT)  [25].  The  majority  of  digital  systems  that  exploit 
the  spectral  content  of  a  signal,  however,  sample  the  signal  and  then  by  processing  it 
appropriately  generate  a  frequency  representation  that  is  also  discrete.  This  trans¬ 
form  is  called  the  Discrete  Fourier  Transform  (DFT)  and  is  dehned  in  Eq  (2.24)  and 
Eq  (2.25).  Recall  that  x^  denotes  the  element  of  the  vector  x  [25]. 


Ns-l 

Xfc  =  ^  xz 

(2.24) 

1=0 

,  Ns-l 

fc=0 

(2.25) 
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where  Ng  is  the  total  number  of  samples.  In  this  case,  however,  the  elements  of  x 
are  samples  of  the  continuous  time  function  x{t)  taken  at  a  sampling  interval  to,  as 
such  Xfc  =  to  x{ktQ).  The  Shannon  Sampling  Theorem  describes  the  circumstances  in 
which  X  completely  represents  x(t).  Sampling  a  continuous  time  signal,  at  to,  results 
in  a  periodic  frequency  spectrum.  The  DTFT  of  the  sampled  signal  is  [25] 

OO  OO 

tox{kto)e-^^^f^^^  =  Y,  X{f  +  k^o).  (2.26) 

k  =  —  00  k=  —  QC 

Therefore,  by  sampling  at  an  interval  to  the  spectrum  of  the  signal  becomes  periodic 
at  $0  =  1/^0-  Effective  reconstruction  of  the  original  signal  is  attained  if  the  spectral 
copies  created  by  sampling  are  removed.  Therefore,  two  things  are  required  to  be 
able  to  reconstruct  x{t),  hrstly  the  spectral  copies  must  not  interfere  with  each  other, 
and  secondly  each  of  the  copies  for  k  jtz  Q  must  be  removed.  The  hrst  of  these  is 
achieved  by  ensuring  that  x(t)  is  band-limited,  that  is  x(t)  G  S 3(^2 f max),  and  to 
must  be  chosen  such  that  $0  =  l/^o  >  2fmax,  be.,  to  avoid  aliasing  the  Nyquist  rate 
must  be  achieved  [25:129].  Secondly,  removal  of  the  spectral  copies  of  the  sampled 
signal  is  accomplished  by  passing  only  those  frequencies  |/|  <  fmax  <  *ho/2.  When 
*^0  =  2 fmax  an  ideal  low-pass  hlter  is  required,  dehned  in  Eq  (2.27),  with  impulse 
response  h(t)  described  as 


where 


H{f)  =  rect(//$o) 


=  — sinc(<I>of), 

to 


rect(0 


1  e<l/2 

0  else 


(2.27) 


(2.28) 


(2.29) 


2-9 


and 


sinc(0 


(2.30) 


sin(7rO 


Therefore  reconstruction  of  x{t)  is  possible  via  the  convolution  of  the  sampled 
signal  with  the  impulse  response  function  h{t): 


x{t)  =  tf)  xikto) -k  hit) 

OO 

=  tQx{kto)  hit  —  kto) 

k=—oo 

OO 

=  a;(A;to)  sinc(<ho[^  —  ^^o])-  (2-31) 

k=—oo 

In  reality,  h{t)  is  not  a  time-limited  signal,  therefore  this  type  of  reconstruction  is 
not  realizable.  In  practice  to  is  chosen  so  that  $0  3>  ‘If max  enabling  the  use  of  cheap 
hlters  with  smooth  pass-band  characteristics,  he.,  H{f)  ^  Yect{f /2fmax)- 

The  transition  to  sampled  signals  from  continuous  time  is  therefore  permissable 
provided  the  Nyquist  sample  rate  is  achieved.  In  numerically  approximating  the 
Fourier  Transform  of  a  sampled  signal,  however,  further  constraints  are  encountered. 
Beginning  by  dehning  the  rectangle-rule  approximation  to  the  Fourier  transform 
consider  signals  sampled  in  both  time,  sample  interval  to,  and  frequency,  sample 
interval  /q. 

/OO 

■OO 

^  X{lfo)  e^2-To*j^  ^  (2.32) 

/=o 

This  approximation  leads  to  the  mapping  of  discrete-frequency  signals  to  discrete¬ 
time  signals.  Having  Ng  samples  in  frequency  leads  to  a  transform  that  calculates 
the  corresponding  iVs-sample  time  signal,  tox{kto).  From  Eq  (2.26)  $0^0  =  1;  and 
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since  there  are  Ng  samples  in  frequency  at  /o,  $0  =  Ngfo  to  fo  =  l/Ng.  The 
sampled  time  signal  is  therefore  represented  by 

Ns-l 

toxikto)  =  toY.X{lfo)e^‘^-^^/^^fo 

1=0 

.  Ns-l 

=  ^  X(//o)  (2.33) 

/=o 

which  is  the  DFT  as  dehned  in  Eq  (2.24).  The  choice  of  Ng  and  the  associated 
sample  intervals  dehnes  the  resolution  of  the  transform  and  extent  of  aliasing.  As 
shown  earlier,  to  prevent  aliasing  1/to  ^  ‘^fmax-  The  sample  interval  and  rate,  to 
and  /o,  are  related  via  to  fo  =  ^/Ng.  Since  to  must  be  chosen  to  prevent  aliasing, 
the  spectral  sample  interval,  fo,  can  only  be  shortened  by  increasing  the  number  of 
samples  Ng,  a  technique  known  as  zero  padding. 

2.3  Pulse  Doppler  Radar  Fundamentals 

2.3.1  Radar  and  Communieations  Model.  This  section  dehnes  the  radar 
model  that  is  assumed  in  this  thesis.  In  dehning  the  model  composition  an  arbitrary 
radar  signal  is  considered  and  the  environmental  range  map  that  it  generates  is  ana¬ 
lyzed.  Radar  and  communications  systems  are  constrained  to  transmit  real  signals. 
However,  the  transmission  of  complex  waveforms  may  be  considered  possible  through 
quadrature  receiver  design  [6].  An  arbitrary  signal  Xmit)  =  a(t)  cos[27r/ct  -|-  fit)], 
has  envelope  a{t),  phase  modulation  (/(t)  and  is  modulated  with  carrier  frequency 
/c  Hz.  When  this  signal  is  processed,  as  illustrated  in  Figure  2.1,  via  mixing  with  an 
in-phase  signal  2  cos{2nfct)  then  low-pass  hltered  the  signal  a(t)  cos[0(t)]  is  obtained. 
Likewise,  when  the  signal  is  quadrature  mixed  with  — 2  sin(27r/ct),  then  hltered,  the 
result  is  a(t)  sin[0(t)]. 

The  output,  or  low-pass  signal  corresponding  to  Xmit),  may  now  be  described 
in  terms  of  an  ordered  pair  {a(t)  cos[0(t)],  a(t)  sin[0(t)]},  and  represented  in  complex 
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a(t)  cos[(;/)(t)] 


a{t)  sin[(;/)(t)] 


Figure  2.1.  In-phase  and  Quadrature  signal  processing.  The  com¬ 
ponent  labelled  —90°  effects  a  phase  shift  on  the  passed 
signal. 


notation  as  shown  in  Eq  (2.34)  [6]: 


x{t)  =  LFF[xm{t)2cos{27rfct)]+ jLFF[-Xm{t)2sm{27rfct)] 

=  a{t)  cos[0(t)]  -|-  j  a{t)  sin[0(t)] 

=  (2.34) 


Accordingly,  this  notation  is  employed  to  represent  the  signal’s,  Xm(t),  low-pass  in- 
phase  and  quadrature  (IQ)  components,  denoted  the  complex  low-pass  video  signal 
[6:131].  This  means  of  modulation  is  a  form  of  quadrature  multiplexing  and  is  known 
to  suffer  performance  degradation  from  a  loss  of  phase  lock.  The  degradation  occurs 
in  the  form  of  undesirable  cross-talk  in  the  output  IQ  channels.  Radar  systems  that 
employ  complex  low-pass  video  signals,  therefore  require  coherent  receiver  and  trans¬ 
mitter  coupling.  Assuming  that  the  reference  radar  system  is  capable  of  coherent 
processing  and  immobile,  a  single  sweep  of  the  range  channel  is  analyzed. 

The  electromagnetic  energy  returned  to  a  radar  at  any  given  time  is  deter¬ 
mined  by  the  spectral  reflectivity  of  a  spherical  range  shell  about  the  emitting  radar 
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antenna  [6:136].  The  reflectivity  of  the  environment  is  modelled  with  a  function, 
u{s),  representing  the  radial  reflectivity  density  defined  by 


u{s)  =  UQ{s)e^‘^^^\ 


(2.35) 


Subsequently,  the  incremental  radar  return  in  the  range  channel  due  to  signal  Xm(t), 
illuminating  a  shell  at  range  s  and  thickness  ds,  is  represented  by: 

Uo{s)ds  a(t  — 


2s 


cos 


27r/c  t  - 


2s 


+ 


2s  .  ,  , 

t - +  ip[s) 

c 


(2.36) 


where  c  is  the  speed  of  light,  and  2s /c  represents  the  round-trip  time  delay  between 
transmission  and  reception  of  Xm{t)  from  a  reflector  at  range  s.  The  coherent  super¬ 
position  of  these  returns  describes  the  radar  return  in  the  range  channel,  calculated 
here  as  [6]: 


Tmit) 


cos 


27r/c 


+  0 


+  (^(s) 


ds. 


(2.37) 


The  previous  analysis  of  IQ  demodulation  of  Xmit)  lead  to  the  compact  complex 
low-pass  video  notation  x{t).  Therefore,  when  r-mit)  is  coherently  IQ  demodulated 
Eq  (2.37)  becomes 


r{t) 


(2.38) 


With  a  simple  change  of  variables,  the  delay  term  is  redefined  as  2s /c  =  r  and  the 
form  of  a  convolution  becomes  more  readily  apparent.  Finally,  Eq  (2.38)  becomes 


r{t) 


_T)eJ‘Pd-r)_dT 


Um(t)  -k  x(t). 


(2.39) 
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The  radar  range  channel  can  therefore  be  represented  by  the  convolution  of  the  mod¬ 
ulated  radial  reflectivity  density  with  the  complex  low-pass  video  of  the  transmitted 
signal.  This  model  is  illustrated  in  Figure  2.2.  To  more  accurately  capture  the  effects 
of  environmental  interference  the  range  cell  is  combined  with  additive  noise.  The 
hnal  component  of  the  model  is  the  receiver  hlter. 


Reflectivity 
density  u{t) 


n(t) 

1 

x(t) 

»■(*) 

h{t) 

X{f) 

HU) 

output 


Propagation 
phase  shift 


x{t)  Complex  modulation  of  the  transmitted  pulse 
h(t)  Impulse  response  of  receiver 


Figure  2.2.  Radar  model  [6:138]  relating  spectral  reflectivity  density 
u{s),  the  complex  modulation  of  the  transmitted  pulse 
x(t)  and  radar  return. 


Note  that  the  true  reflectivity  density  is  effected  by  the  radar  waveform  se¬ 
lection,  modulation  frequency,  and  the  receiver  hlter  design.  These  aspects  of  the 
model  are  discussed  in  the  following  section. 

2.3.2  Bandpass  Signals  and  Systems.  This  section  revisits  some  of  the 
characteristics  of  bandpass  signals  and  systems  that  are  key  to  the  development  of 
this  thesis.  Due  to  platform  size  limitations  and  antenna  gain  characteristics^  radar 
signals  are  typically  modulated  at  high  frequencies.  Also,  due  to  power  constraints 
and  spectrum  utilization  these  signals  are  bandlimited.  Subsequently,  the  trans¬ 
mitted  radar  signal  is  described  as  narrowband,  meaning  the  carrier  frequency  is 
much  greater  than  the  signal  bandwidth,  i.e.,  {x(t)  :  x(t)  G  with  /q  S>  $. 

Primarily,  this  permits  two  things  that  are  discussed  in  this  section:  hrstly,  the  sim- 
plihed  use  of  analytic  notation,  and  secondly  optimization  of  receiver  hlter  design  to 
eliminate  noise. 

^For  example,  the  gain  of  a  parabolic  antenna  is  inversely  proportional  to  the  square  of  the 
wavelength  of  the  transmitted  signal,  recall  G  «  [35] . 
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The  analytic  representation  of  a  signal  'ip(t)  =  x(t)  +  jx(t),  from  Eq  (2.19), 
describes  a  signal  having  only  a  positive-sided  spectrum  and  real  component  x{t), 
however  the  calculation  of  the  Hilbert  transform  can  generate  unwieldy  functions. 
Notably,  if  x{t)  is  modulated  with  a  complex  carrier,  its  entire  spectrum  is  shifted 
from  baseband  to  fc,  so 

Xa{t)  = 

=  a{t)  cos[27r/ct  -|-  4>{t)]  +  j  sin[27r/ct  -|-  0(t)] .  (2.40) 

If  the  shifted  spectrum  has  only  positive  frequency  content  then  it  approximately 
meets  the  criteria  for  the  analytic  form  for  Xm{t)-  For  this  to  be  analytic  in  the  strict 
sense,  however,  the  imaginary  part  of  Xa{t)  must  be  the  Hilbert  transform  of  its 
real  part.  A  detailed  analysis  of  the  Hilbert  transform  of  Xm{t)  shows  that  provided 
=  a(t)  cos[0(t)]  and  A{a;(t)}  =  a(t)  sin[(;/)(t)]  are  bandlimited  below  fc,  then 
Xa(t),  from  Eq  (2.40),  is  the  correct  analytic  signal  representation  of  x^it)  [4]. 

Figure  2.3  shows  the  effect  of  modulating  a  signal  on  a  complex  carrier.  Fig¬ 
ure  2.3(b)  illustrates  the  spectrum  modulated  insufficiently,  subsequently  containing 
negative  frequencies  and  therefore  failing  as  an  analytic  approximation.  In  Figure 
2.3(c)  the  modulation  frequency  is  chosen  such  that  it  is  greater  than  the  half¬ 
bandwidth  of  the  signal.  Radar  signals  are  typically  considered  narrowband,  and  are 
therefore  able  to  utilize  Eq  (2.40)  as  a  valid  analytic  representation  of  Xm(t). 

The  receiver  design  for  narrowband  signals  is  next  discussed  by  considering  the 
detection  of  a  single  target  in  stationary  white  Gaussian  noise.  Although  optimality 
in  this  scenario  depends  on  noise  parameters  and  costs  associated  with  detection 
threshold  selection,  the  following  investigation  considers  receiver  optimization  in 
terms  of  maximizing  only  the  Signal-to-Noise  Ratio  {SNfZ).  Figure  2.2  provides  the 
model  in  this  case,  by  considering  a  signal  in  additive  noise  y  =  r  +  n.  The  system 
output,  therefore,  is  z  =  y'kh  =  r-kh  +  n'kh.  The  SAfTZ  under  consideration  is 
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Figure  2.3.  Spectral  modulation  of  a  bandlimited  signal.  A  bandlim- 
ited  signal,  (a),  can  be  represented  using  analytic  signal 
notion  provided  its  modulated  spectrum,  (b),  exists  in 
only  one  hemisphere,  (c). 


dehned  as  the  ratio  of  peak  signal  power  to  average  noise  power: 


1“  S„(/)|i/(/)|V/ 


(2.41) 


where  S'„(/)  is  the  power  spectral  density  of  n{t).  By  the  Schwartz  inequality  the 
numerator  can  be  bounded  as 


< 


\R{f)\^df  /  \H{f)e^^-f-\^df 


’  — oo 
poo 


’  —oo 
poo 


\R{f)\^df  /  \H{f)\^df, 


(2.42) 
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assuming  n{t)  is  white  Gaussian  noise  with  power  spectral  density  S'„(/)  =  rj/2. 
Combining  this  with  Eq  (2.41)  and  Eq  (2.42)  the  SMTZ  becomes 


ssfn  < 


-  /  \RU)?df 

V  J-oo 

28 

V  ' 


(2.43) 


Equality  holds  in  Eq  (2.43),  if  H{f)  =  kR*{f)e~^‘^'"R .  Therefore  the  optimal  filter 
in  white  Gaussian  noise  is  the  matched  filter,  so  named  because  the  spectrum  of 
the  receiver  is  matched  to  that  of  the  noise  free  signal.  Consequently,  the  impulse 
response  function  of  H{f)  is  h(t)  =  kr*{T  —  t),  a  time  delayed  and  flipped  conjugate 
of  the  noise  free  signal. 

A  single  point-target  at  range  R  and  with  radial  velocity  Vr  has  a  two-fold 
effect  on  the  returned  constant  carrier,  /c,  pulse.  The  range  of  the  target  contributes 
to  a  round  trip  delay  of  r  =  2i?/c,  and  the  radial  velocity  results  in  a  Doppler 
shift,  V  =  2vrfc/c.  The  output  of  the  matched  filter  can  then  be  approximated  for  a 
single-target  response  as  a  function  of  Doppler  and  delay,  as  detailed  in  Eq  (2.44): 


x(r,  z/)  =  /  x{t)x*{t  — 


(2.44) 


This  function  was  formulated  by  Ville  [40]  and  introduced  to  radar  by  Woodward  [41] 
and  is  commonly  known  as  the  Ambiguity  Function  [31].  Though  modelling  a  moving 
target  as  a  spectral  shift  is  over  simplified,  this  function  serves  as  a  good  estimate 
in  describing  the  ambiguities  that  exist  for  a  given  transmitted  waveform.  When 
normalized,  this  surface  provides  a  quasi-joint  density  for  describing  the  probability 
of  a  target’s  range  and  velocity.  Appendix  A  provides  ambiguity  surfaces  for  a 
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single  pulse  at  Appendix  A.  1.1,  a  pulse  train  at  Appendix  A.  1.2  and  a  LFM  pulse 
at  Appendix  A.  1.3. 

2.3.3  Radar  Waveform  Design.  The  previous  sections  describes  the  signal 
characteristics  and  radar  model,  as  shown  in  Figure  2.2.  This  naturally  leads  to 
the  topic  of  signal  design.  The  constraints  limiting  the  antenna  size  also  restrict 
the  number  of  antennae.  Typically,  a  single  antenna  is  used  in  the  majority  of 
airborne  platforms,  which  logically  constrains  the  analysis  to  monostatic  systems, 
he.,  collocated  transmitter  and  receiver.  Such  systems  must  use  pulsed  waveforms  to 
provide  a  listening  interval  to  detect  radar  returns.  This  section  explores  desirable 
properties  of  radar  waveforms  and  the  design  trade-offs  associated  with  inherent 
ambiguities  and  resolution. 

2.3.3. 1  Pulse  Doppler  Radar  Ambiguities.  A  radar  is  primarily  em¬ 
ployed  for  determining  information  about  target  presence,  and  associated  parameters 
such  as  size,  location  and  motion.  This  information  is  delivered  in  a  hierarchy  from 
target  detection  to  parameter  estimation  with  the  governing  metric  of  precision  a 
function  of  resolution  [31].  To  improve  a  radar’s  performance  in  target  detection, 
the  SJ\fTZ  is  maximized  using  a  matched  hlter  as  described  in  Section  2.3.2.  For 
parameter  estimation,  a  target’s  range  and  velocity  is  considered  via  analyzing  the 
matched  hlter  output. 

By  dividing  the  radar  range  channel  into  bins,  an  estimate  of  range  is  given  by 
detecting  the  presence  of  a  pulse  return  in  a  given  bin.  If  a  pulse  x(t)  =rect(t/r),  as 
shown  in  Figure  2.4(a),  were  to  be  transmitted  then  the  matched  hlter  output  would 
be  a  triangular  pulse  with  width  2r  as  shown  in  Figure  2.4(b).  Using  the  half-power 
point  as  a  measure  of  pulse  width,  the  range  bins  become  r  wide  and  provide  a  range 
resolution  of  cr.  Target  velocity,  on  the  other  hand,  is  determined  as  the  range  rate 
and  can  be  calculated  from  the  Doppler  shift  of  the  return  signal  [35].  The  original 
spectrum  of  the  pulse  is  a  sinc(-)  waveform,  but  the  spectrum  of  the  matched  hlter 
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Figure  2.4.  Matched  filter  output  of  (a)  in  time,  (b),  and  frequency, 
(c).  The  matched  hltering  of  a  rect-pulse  results  in  a  de¬ 
crease  in  time  resolution  but  an  increase  in  the  frequency 
resolution. 


output  is  a  sinc^(-)  waveform,  as  shown  in  Figure  2.4(c).  Consequently,  the  matched 
hlter  output  is  narrower  in  frequency  than  the  original  signal  providing  better  fre¬ 
quency  resolution,  illustrating  a  clear  relationship  between  detection,  estimation  and 
resolution. 

However,  pulses  are  transmitted  periodically,  therefore  the  transmitted  signal 
is  described  in  Eq  (2.45),  where  T  is  the  Pulse  Repetition  Interval  (PRI)  and  r  is 
the  pulse  width: 

x(t)  =  ^  rect  (- — (2.45) 

n=—oo  ^  ^  • 

This  leads  to  two  signihcant  observations.  Firstly,  if  the  radar  power  is  such  that  the 
maximum  detection  range  is  greater  than  cT/2,  then  the  range  returns  are  ambigu¬ 
ous.  Figure  2.5  illustrates  this  showing  that  if  two  targets,  at  ranges  R1  and  R2,  are 
being  illuminated  by  a  pulsed  radar  then  the  returns  from  targets  in  the  ambiguous 
range  region,  R  >  Ru  =  cT/2,  fold-over  and  generate  misleading  range  indications. 
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Secondly,  the  periodicity  of  the  signal  results  in  a  sampling  of  the  signal  spectra^. 
The  spectrum  of  a  T-periodic  pulse  train  is  simply  that  of  a  single  pulse  sampled 
at  an  interval  1/T.  Therefore,  if  a  pulse  train  is  processed  over  what  is  termed  a 
Coherent  Processing  Interval  (CPI),  then  a  comb-like  spectra  is  generated  resulting 
in  better  Doppler  resolution.  However,  because  the  spectra  is  evenly  spaced  at  1/T, 
ambiguities  in  the  measurement  of  spectral  shifts,  or  Doppler,  are  introduced.  The 
unambiguous  radial  velocity  of  a  target  is  |n|  <  |n„|  =  A/4T.  Therefore,  unambigu¬ 
ous  radar  operation  requires  i?|n|  <  cA/8. 
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Figure  2.5.  Ambiguous  pulsed  radar  returns  [1:1-3].  The  target  at  R2 
is  located  in  the  ambiguous  range  interval,  subsequently 
the  hrst  pulse,  PI,  is  processed  in  the  second  CPI  as  a 
return  from  pulse  P2. 


Obviously,  a  better  estimate  of  target  parameters  is  attained  with  the  removal 
of  all  types  of  ambiguities,  and  increased  resolution  in  both  time  and  frequency. 
Notably,  by  increasing  the  PRI,  the  unambiguous  range  is  increased.  However,  the 
increase  in  PRI  conversely  decreases  unambiguous  velocity.  Clearly,  simultaneous 
removal  of  ambiguities  in  both  range  and  Doppler  cannot  be  achieved  by  solely 
altering  the  PRI.  A  technique,  know  as  Non-linear  Ambiguity  Suppression  (NFS), 
solves  this  problem  to  a  considerable  extent  and  will  be  discussed  in  Section  2.4. 

^Recall  from  a  Fourier  series  analysis  that  if  a  signal  is  T  periodic  then  it  can  be  represented 
by  xrit)  =  Y^^=-oo  where  /o  =  1/T  and  =  foX{nfo).  Therefore,  the  spectrum  of 

XT{t)  becomes  T{xT{t)}  =  J2T=-oo  foX{nfo)6{f  -  nfo)  [39]. 
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2. 3. 3. 2  Pulse  Compression.  Previously,  the  resolution  of  the  radar 
pulse  was  arbitrarily  dehned  as  the  half-power  region  of  a  signal;  a  more  rigorous 
analysis  of  resolution  is  now  presented.  A  metric  is  required  that  operates  on  a  signal 
and  provides  a  real  number  describing  the  resolution  of  that  signal,  a  :  {x}  — ^  TZ. 
Three  such  metrics  are:  total  duration  (ur),  equivalent  rectangle  {ce),  and  root 
mean  squared  (rms)  signal  duration  (a).  Total  duration  is  dehned  in  time  and 
frequency  respectively  as: 

(2.46) 


The  equivalent  rectangle  metric  aims  to  describe  the  width  required  to 
dehne  a  rect(-)  pulse  of  equivalent  height  to  the  signal  at  its  center  such  that  the 
rect(-)  contains  the  same  amount  of  energy  as  the  signal.  This  is  dehned  in  time  and 
frequency  respectively  as: 


o'e 


fZc  " 


(2.47) 


Finally  the  rms  signal  duration  is  dehned.  The  moments  of  statistical  densities 
provide  information  about  spread  and  location  and  are  therefore  used  to  characterize 
random  variables.  Analogously,  these  moments  are  often  used  to  characterize  signal 
precision  and  resolution  performance  [31].  The  hrst  moments  are  used  to  describe  a 
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signal’s  center  of  mass  and  are  defined  in  time  and  frequency  respectively  as: 


/!!  \x{t)\^dt 

(2.48) 

IZof\Xif)M 


The  centered  second  moment,  denoted  here  as  rms  signal  duration,  is  a  measure  of 
spread  and  is  dehned  in  time  and  frequency  respectively  in  Eq  (2.49).  These  metrics 
are  used  most  predominantly  throughout  this  thesis. 


(2.49) 

irjf  - 

JZ\X{fWdf 


By  invoking  Schwartz  inequality  on  the  product  of  the  time  and  frequency  rms 
signal  durations  of  Eq  (2.49),  a  time-bandwidth  product  is  derived  and  bounded  by 
Gabor  [12]  so  that 

aS  >  1,  (2.50) 

where  equality  is  achieved  when  x(t)  =  — [31:55].  As  an  initial  observation, 
this  indicates  that  by  increasing  signal  bandwidth  the  rms  time  duration  reduces  and 
therefore  range  resolution  improves  proportionately,  with  a  LEM  pulse  providing  the 
best  performance. 

This  indicates  that  the  range  resolution  of  a  signal  is,  at  the  very  least,  in¬ 
versely  proportional  to  its  bandwidth,  providing  an  upper  bound.  By  exploiting  this 
property  and  making  the  time-bandwidth  product  as  large  as  the  system  can  sustain, 
the  greatest  resolution  for  a  given  bandwidth  is  achieved.  For  uniform  pulsed  signals, 
a  popular  means  of  doing  this  is  through  LEM.  Such  pulses  are  generated  by  either 
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directly  modulating  a  constant  frequency  signal  with  a  variable-frequency  oscillator, 
or  passing  a  constant  frequency  signal  through  a  dispersive  hlter  [17].  These  signals 
are  known  as  either  LFM  or  quadratic  phase  signals  having  the  form  described  in 
Eq  (2.51): 

a^LFM(t)  =  rect^-^  where  p  =  <h/r.  (2-51) 

The  matched  hlter  output  of  this  signal  is 

,  ,  f  t\  sin7r[utr  — 

(2.52) 


characterized  by  quadratic  phase  and  with  its  maximum  at  t  =  0. 

An  understanding  of  how  bandwidth  affects  the  resolution  of  this  function  is 
obtained  by  considering  the  null-null  width,  a,  of  the  mainlobe: 


A/)Ur^-4 

2v/p 


(2.53) 


When  the  pulse  width  r  is  hxed,  the  bandwidth  $  increases  linearly  with  p.  Therefore 
as  fi  is  increased,  and  thus  the  bandwidth  of  the  signal  increased,  the  numerator 
approaches  zero,  he.,  a/ /ir^  -|-  4  ^  a/ /ir^  —  4,  and  the  denominator  increases,  causing 
a  to  become  very  small.  LFM  is  a  means  of  improving  the  resolution  of  a  matched 
hlter  output  by  increasing  the  spectral  content  and  bandwidth  of  a  hxed  duration 
single  pulse. 

In  general,  a  waveform  that  has  a  time-bandwidth  product  much  greater  than 
unity  is  said  to  exhibit  properties  of  pulse  compression  [31].  There  are  many  wave¬ 
forms  that  have  been  designed,  aside  from  LFM,  exhibiting  these  properties  such  as 
Barker  codes  [35],  polyphase  codes  [21]  [22],  Costas  codes  [8],  and  maximal  length 
pseudo-random  noise  sequences  [23],  such  as  Gold  codes  [7]. 
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Ambiguity  Suppression 

The  ambiguities  detailed  in  Section  2. 3. 3.1  are  inherent  for  any  pulsed  radar 
system.  This  section  discusses  approaches  to  mitigate  these  ambiguities,  discussing 
a  multiple  PRF  approach  and  a  pulse  discrimination  process  known  as  Non-linear 
Ambiguity  Suppression  (NLS)  [28]. 

2.4-1  Pulse  Repetition  Frequeney  Approaches.  Ambiguities  can  be  miti¬ 
gated  to  some  extent  through  the  choice  of  the  PRI,  however,  as  this  section  shows, 
these  tend  to  be  application  specihc  leading  to  a  search  for  other  alternatives.  Sec¬ 
tion  2.3.3. 1  discusses  the  origins  of  range  and  velocity  ambiguities  illustrating  that 
the  selection  of  the  PRI  T  jointly  effects  velocity  ambiguity  and  range  ambiguity, 
motivating  an  analysis  of  PRI  selection. 

Frequency 
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Figure  2.6.  Pulse  train  spectra  and  PRI.  An  inhnite  length  pulse 
train  with  PRI  T,  (b),  will  have  a  comb-like  spectra  with 
harmonics  separated  by  1/T,  (a). 


Typically,  pulse  doppler  radars  operate  at  one  of  three  PRI  rates:  High, 
Medium,  or  Low.  The  PRI  rate  is  also  known  as  the  Pulse  Repetition  Frequency 
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(PRF)  and  is  defined  as  fp  =  1/T.  High  PRF  radars  have  greater  spectral  sepa¬ 
ration  and  therefore  are  better  suited  for  avoiding  velocity  ambiguities,  but  due  to 
the  associated  small  PRI  suffer  greatly  from  range  ambiguities.  Low  PRF  radars,  on 
the  other  hand,  have  an  unambiguous  range  near  the  maximum  detection  range  of 
the  radar  system,  but  have  a  densely  packed  spectra  exacerbating  velocity  ambigui¬ 
ties.  Figure  2.6(a)  illustrates  the  spectrum  resulting  from  an  inhnite  train  of  pulses, 
shown  in  Figure  2.6(b).  The  aforementioned  accordion  effect  that  PRI  changes  have 
on  the  spectral  separation  is  clearly  represented  in  this  hgure. 

By  utilizing  multiple  PRF’s,  these  ambiguities  can  be  signihcantly  reduced. 
This  technique  may  be  conceptualized  as  the  transmission  of  two  or  more  pulsed 
waveforms  from  the  same  antenna,  with  each  waveform  having  a  unique  PRF.  The 
PRF’s  are  typically  related  via 


PRF.  =  (2.54) 

where  PRF.  is  the  PRF,  a.  is  an  integer,  and  is  a  common  sub-multiple  of 
a  prescribed  PRI  [1:1-18]  [24:273].  A  radar  return  from  a  target  has  a  true  range 
Rt  described  by  Rt  =  yR.  -|-  Rj,  where  Ri  is  the  indicated  range  and  yR.  is  an 
integer  multiple  of  the  unambiguous  range.  The  unknowns  are  y  and  Rt;  therefore, 
by  changing  the  PRF  (and  thereby  R.  and  R*,)  a  second  simultaneous  equation 
is  generated  and  a  solution  for  the  true  range  is  found.  Alternatively,  the  true 
range  can  be  determined  by  identifying  where  the  radar  returns  for  the  different 
PRFs  coincide.  However,  the  multiple  PRF  returns  from  this  technique  cannot  be 
processed  coherently,  resulting  in  a  less  energy  efficient  implementation.  Also,  an 
environment  that  has  multiple  targets  is  likely  to  generate  further  ambiguous  returns 
from  ghosting  [1:1-18]. 

2.4-2  Non-linear  Suppression  System  Model.  Non-linear  Ambiguity  Sup¬ 
pression  (NLS)  combines  pulse  compression  and  pulse  discrimination  techniques  to 
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deliver  a  means  of  improving  radar  ambiguities  in  pulsed  radar.  The  concept  was 
introduced  by  Palermo  in  1962  [28]  utilizing  a  conjugate  LFM  coded  pulse  pair. 
The  fundamental  building  block  of  the  NLS  technique  is  the  suppression  kernel 
shown  in  Figure  2.7.  This  unit  operates  on  the  radar  return  to  suppress  a  particular 
pulse  response  by  hrst  pre-hltering  the  entire  signal  to  discriminate  that  pulse  in 
some  way  from  all  other  pulses.  Once  the  presence  of  that  pulse  is  emphasized  in 
the  radar  return,  the  pulse  is  suppressed  by  a  suppression  operator,  which  ideally 
has  no  effect  on  all  other  pulses.  After  this  suppression  operation,  the  radar  return 
is  passed  through  a  post-hltering  process  to  remove  the  effect  of  pre-hltering.  The 
pre-hlter  hk{t)  and  post-hlter  gk{t)  condition  the  waveform  to  maximize  the  effect 
of  the  suppression  operation  on  the  pulse  without  impacting  the  other  pulses. 
Therefore,  the  suppression  operation  determines  the  hlter  and  pulse  design  criteria. 


Afc 


output 


Figure  2.7.  General  NLS  kernel  [1:2-11]. 


Two  examples  of  thresholding  suppression  operators  are  illustrated  in  Fig¬ 
ure  2.8(a)  and  (b).  Although  various  suppression  operators  may  be  devised,  this 
thesis  focuses  only  on  the  use  of  thresholding  suppression  operators  as  these  are  the 
most  common.  The  implications  that  thresholding  has  on  formulation  of  a  pulse  de¬ 
sign  criteria  is  two-fold.  To  remove  energy  from  a  signal,  the  focused  pulse  response 
must  have  an  amplitude  greater  than  the  threshold,  consequently  the  pre-hltering 
process  should  project  as  much  of  the  undesirable  pulse  energy  above  the  threshold 
as  possible  whilst  spreading  all  other  pulses  below  it.  When  only  thresholding  sup¬ 
pression  operators  are  considered,  the  suppression  is  dependent  on  the  ability  of  the 
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Figure  2.8.  Threshold  suppression  operators  (r^).  Hole-punch  (a), 
and  saturation  limiter  (b),  with  a  =  0.51^^  [1:1-23]. 

pre-hlter  to  compress  the  intended  pulse  whilst  spreading,  or  dispersing,  all  other 
pulses. 

A  common  means  of  pulse  compression  exists  through  the  matched  hltering  of 
certain  waveforms.  To  utilize  the  NTS  technique,  a  series  of  pulses  must  be  designed 
in  snch  a  way  that  as  a  resnlt  of  matched  hltering  for  the  pnlse  (causing  it  to 
become  compressed)  each  of  the  remaining  pulses  are  dispersed.  This  concept  is 
known  as  mntnal  dispersion  and  is  ntilized  in  spread-spectrnm  applications  snch  as 
direct  sequence  Code  Division  Mnltiple  Access  (CDMA)  [30].  When  thresholding 
NTS  is  employed,  the  components  of  the  snppression  kernel,  in  Fignre  2.7,  become: 
matched  hlter,  threshold  operator,  and  conjngate  rehection  hlter. 

A  NTS  system  with  iV  pnlses  ntilizes  as  many  channels  to  disinter  each  pulse 
from  the  radar  return.  The  channel  contains  iV— 1  snppression  kernels,  A^.  V/c  7^  n, 
aimed  at  suppressing  all  but  the  pulse  from  the  pulse  train  as  detailed  in  Fig¬ 
nre  2.9.  The  hnal  stage  of  each  channel  is  a  matched  hlter  operation  for  the  songht- 
after  pulse.  In  other  words,  the  channel  extracts  the  pnlse  from  a  radar 
retnrn  whilst  snppressing  all  others.  As  a  resnlt,  the  pnlse  repetition  interval  ehec- 
tively  increases  iV-fold.  This  extends  the  ehective  nnambignous  range  to  NcT/2  and 
considerably  alleviates  the  range  ambignities  discussed  in  Section  2.3.3. 
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Figure  2.9.  NFS  system  illustrating  N  individual  channels  each  com¬ 
prising  TV  — 1  suppression  operations  [1:2-15]. 

To  illustrate  this  approach  consider  a  pulse  family  comprised  of  two  conjugate 
LFM  waveforms  /i(t)  and  /2(t)  for  |t|  <  1  and  p  =  5.  The 

radar  range  channel  is  arbitrarily  constructed  over  a  CPI  to  contain  these  pulses, 
r(t)  =  flit  —  ti)  +  f2{t  —  t2)  +  n{t),  delayed  in  time  and  affected  by  Gaussian  white 
noise,  as  detailed  in  Figure  2.10(a).  Notice  that  the  pulse  returns  overlap;  this 
phenomena  is  known  as  self-clutter  due  to  the  interference  each  pulse  contributes 
to  the  other.  In  attempting  to  disinter  /2(t)  from  this  signal  it  is  passed  through 
a  matched  hlter,  F^if),  resulting  in  the  plot  in  Figure  2.10(b).  On  the  other  hand 
suppressing  fi{t)  with  a  thresholding  non-linearity  removes  a  good  deal  of  its  energy, 
and  therefore  interference  contribution,  as  shown  in  2.10(c).  Because  of  the  dispersive 
effect  that  matched  hltering  for  /2(t)  has  on  fi(t),  both  approaches  provide  good 
results,  however  the  contribution  to  self-clutter  is  significantly  reduced  due  to  the 
suppression  operation  as  shown  in  Figures  2.10(b)  and  (c). 

2. 5  Summary 

This  chapter  introduced  the  background  theory  used  throughout  the  remain¬ 
der  of  this  thesis,  and  in  so  doing  presents  a  cross-section  of  the  literature  pertinent 
to  this  held.  Section  2.2.1  introduced  the  concept  of  realizability  and  some  signal 
theory  basics  and  metrics.  This  analysis  naturally  lead  to  a  discussion  on  integral 
transforms,  the  most  prevalent  being  the  Fourier  and  Hilbert  transform,  and  to  sub- 
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Figure  2.10.  NFS  processing  of  conjugate  LFM  pulses.  A  radar  re¬ 
turn  comprising  a  pair  of  conjugate  LFM  pulse  responses 
is  simulated,  (a).  The  return  is  matched  hltered  for 
pulse  2  (right),  compressing  it  and  dispersing  pulse  1 
(left),  (b).  Pulse  1  is  suppressed  prior  to  matched  hlter- 
ing  for  pulse  2,  (c). 
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sequent  consideration  given  to  numerical  approximations.  Section  2.3  introduced 
radar  signals  as  a  subset  of  realizable  signals,  narrowband  in  nature,  permitting 
compact  analytic  notation.  This  section  also  provided  the  model  for  the  radar  sys¬ 
tem.  Section  2.3  concluded  with  a  discussion  on  radar  waveform  design  describing 
ambiguities  inherent  to  pulsed  radar  systems  and  dehning  metrics  for  measuring  res¬ 
olution.  The  core  topic  of  this  chapter  was  introduced  in  Section  2.4.2  where  the 
NTS  approach  was  explained.  The  NTS  model  was  provided  and  NLS’s  superior 
ambiguity  suppression  performance  was  demonstrated. 
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III.  NLS  Symbol  Design  and  Modelling 

3. 1  Introduction 

The  effectiveness  of  a  thresholding  NLS  system  hinges  on  the  selection  of  sym¬ 
bols  that  exhibit  good  compression  in  autocorrelation  and  dispersion  in  crosscorrela¬ 
tion.  Section  2.4.2  provided  an  example  using  conjugate  linear  FM  symbols,  unfor¬ 
tunately  the  use  of  these  symbols  is  limited  to  only  one  pair,  therefore  to  expand  the 
symbol  set  a  more  detailed  analysis  of  NLS  is  required.  This  chapter  describes  an 
approach  conceived  by  Brown  [5]  that  is  optimal  in  rms  time  duration  for  generating 
a  family  of  mutually  dispersive  symbols  of  any  size.  Section  3.2  provides  a  brief 
overview  of  the  NLS  concept.  Section  3.3  discusses  the  scope  of  the  research,  relat¬ 
ing  the  approach  and  associated  synthesis  constraints.  Section  3.4  relates  the  radar 
model  to  NLS  by  dehning  analytically,  as  in  [5],  the  problem  of  designing  a  symbol 
set  suited  to  NLS,  and  illustrates  the  Brown  NLS  system.  Section  3.5  derives  Brown’s 
theorem  with  the  optimal  envelope  and  phase-rate  functions,  introduces  the  Hermit 
Problem  and  dispersive  gain,  and  provides  an  expression  for  theoretically  optimal 
performance.  Section  3.6  dehnes  notation  for  biphase  codes  relating  the  suitability 
of  these  codes  to  NLS.  Finally,  Section  3.7  discusses  the  metrics  used  to  characterize 
and  contrast  Brown  symbols,  providing  an  overview  of  the  approach  to  numerical 
approximation. 

3.2  Non-linear  Ambiguity  Suppression  Overview 

Radar  systems  are  employed  to  detect  targets  and  to  facilitate  parameter  esti¬ 
mation.  Within  this  construct,  the  pulse-like  nature  of  pulse  Doppler  radar  systems 
leads  to  radar  returns  containing  range  and  velocity  ambiguities  degrading  the  over¬ 
all  system  performance.  The  ambiguities  may  be  eliminated  should  the  attitude  of 
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the  target  in  question  be  such  that  the  range  velocity  product  is  [35] 

R\v\  <  ^cA,  (3.1) 

where  c  is  the  speed  of  light  and  A  is  the  wavelength  of  the  transmitted  radar  sig¬ 
nal.  Although  a  prohle  of  the  target’s  range  and  velocity  characteristics  may  be 
constructed,  typically  they  cannot  be  controlled,  requiring  an  alternative  means  of 
eliminating  ambiguities. 

Non-linear  Ambiguity  Suppression  (NLS)  is  a  means  of  transmitting  a  pulse- 
train  comprised  of  N  pulses,  each  of  which  exhibit  specihc  characteristics  derived 
from  the  choice  of  suppression  operator.  When  considering  thresholding  suppression 
operators  the  pre-hltering  the  pulse  train  must  provide  a  two-fold  affect.  The  hlter 
must  compress  the  intended  pulse  and  spread,  or  disperse,  all  remaining  pulses.  This 
is  the  key  design  criteria  for  NLS  symbol  design. 

By  concatenating  a  chain  of  suppression  kernels,  each  suppressing  a  different 
pulse  as  detailed  in  Figure  2.9,  a  series  of  non-linearities  are  effected  on  the  radar 
return  reducing  the  self-clutter.  Figure  3.1  shows  a  time  slice  series  for  a  LFM  signal 
passed  through  a  single  suppression  kernel.  Figure  3.1(a)  is  the  original  linear  FM 
pulse.  This  signal  is  passed  through  a  matched  hlter;  the  response  is  shown  as  a 
dashed  line  in  Figure  3.1(b),  while  the  solid  line  indicates  the  resultant  signal  after 
hole-punch  suppression.  Finally,  the  suppressed  signal  is  passed  through  a  conjugate 
rehection  hlter  where  the  output,  displayed  in  Figure  3.1(c),  exhibits  a  signihcant 
reduction  in  signal  energy. 

The  suppression  kernel,  A^,  considered  in  this  thesis  takes  the  form  illustrated 
in  Figure  3.2  [1:2-11].  The  pre-hltering  is  performed  by  a  matched  hlter,  compressing 
the  pulse  and  spreading  all  others.  The  suppression  operation  is  accomplished 
with  a  thresholding  non-linearity,  and  hnally  the  post-hltering  is  performed  by  a 
conjugate  rehection  hlter.  Using  this  kernel  constrains  the  search  for  optimality. 
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Figure  3.1.  Thresholding  example  with  LFM  (a)  original  signal,  (b)  thresh¬ 
olding  comparison,  and  (c)  suppressed  signal. 
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Figure  3.2.  Threshold  NFS  kernel  [1:2-11]. 
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Brown  has  devised  a  theorem  that  poses  an  optimal  solution  to  pulse  design 
for  this  conhguration  of  NLS  system  [5] .  Brown  [5]  describes,  in  terms  of  rms  dura¬ 
tion,  the  optimal  correlation  waveform  introducing  an  elegant  notation  simplihcation. 
Brown  considers  the  design  of  phase  functions  such  that  a  set  of  codes  may  be  devised 
with  optimal  mutually  dispersive  properties.  This  step  leads  to  a  problem,  coined 
by  Brown  as  the  Hermit  Problem,  in  which  an  iV-dimensional  solution  is  determined 
by  choosing  N  vectors  in  iV-space  that  are  each  maximally  equidistant  from  each 
other. 

3.3  Scope  of  Res eareh 

As  stated  in  Section  1.2,  the  goals  of  this  thesis  are  to: 

1.  Explain  the  formulation  of  Brown’s  theorem  [5]  whilst  introducing  a  simplihed 
linear  algebraic  notation; 

2.  Describe  a  systematic  process  for  designing  Brown  symbols  through  an  analysis 
of  Brown’s  work  and  by  theoretical  extension; 

3.  Provide  a  theoretical  analysis  for  the  digitization  of  Brown  symbols; 

4.  Evaluate  performance  of  Brown  symbols,  comprising  dispersive  gain,  through 
simulation  and  numerical  approximation;  and 

5.  Compare  the  performance  of  various  Brown  symbols  with  discrete  codes  through 
simulation. 

The  simulations  conducted  in  generating  and  characterizing  Brown  symbols  utilized 
Matlab®.  This  approach  constrains  the  simulation  to  purely  digital  analysis,  provid¬ 
ing  the  theoretical  background  required  to  analyze  Brown  symbols  for  digitization 
and  obtaining  results  using  numerical  approximation  techniques.  Analytic  deriva¬ 
tions  of  closed-form  spectral  dehnitions  of  Brown  symbols  are  determined  and  syn¬ 
thesized  as  frequency  sampled  signals,  the  respective  sampled  time  functions  are  de- 
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termined  with  the  rectangle-rule  approximation  technique.  Finally,  the  time-sampled 
correlation  functions  are  calculated  and  analyzed. 

3.4  Brown  Theorem  Motivation 

Pulsed  radar  systems  are  challenged  to  divine  the  nature  of  their  environment 
by  processing  radar  returns  amidst  self-clutter  and  ambiguities.  Section  2.3.3. 1  high¬ 
lighted  the  causes  of  ambiguities,  concluding  that  it  is  not  possible  to  entirely  remove 
ambiguities  using  only  the  PRF  as  a  control  variable.  Alternatively,  NFS  provides  a 
powerful  technique  to  overcome  these  challenges,  however  it  poses  different  problems 
such  as  symbol  selection,  and  difficulties  in  optimal  symbol  design. 

The  discovery  of  symbol  families  that  are  well-suited  to  NFS  is  not  by  any 
means  a  recent  achievement.  CDMA  schemes  are  an  example  of  the  use  of  discrete 
codes  that  exhibit  the  very  properties  desired  in  NFS  [2].  Instead,  the  search  for 
optimal  symbols  has  historically  been  relegated  to  the  search  for  optimal  discrete 
codes;  an  avenue  of  research  that  has  encompassed  a  vast  array  of  techniques  from 
heuristic  approaches  to  neural  networks  but  that  has  not  fully  unveiled  a  solution 
to  optimality  in  NFS.  By  considering  the  rms  time  duration  of  the  correlation 
functions.  Brown  proposes  an  elegant  theorem  for  the  design  of  optimally  mutual 
dispersive  codes  leading  to  an  expression  for  the  optimal  performance  of  such  codes. 

This  section  presents  Brown’s  theorem  by  hrst  constructing  the  environmental 
model  and  presenting  relevant  assumptions  as  necessary.  The  problem  dehnition  is 
then  stated  and  the  theorem  is  explained.  Finally,  the  resultant  theoretical  optimal 
performance  in  rms  time  duration  is  presented. 

3.4.1  Environment  Construetion.  Pulse  Doppler  radars  are  employed  on 
a  vast  array  of  sensor  platforms.  Of  course,  platform  specihcs  determine  the  partic¬ 
ular  performance  prohle  and  associated  catalysts  to  system  degradation.  However, 
the  problem  of  self-clutter  and  range  and  velocity  ambiguities  appears  somewhat 
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ubiquitous.  This  derivation  utilizes  analytic  notation  and  the  radar  model  described 
in  Section  2.3  and  illustrated  in  Figure  2.2.  Accordingly,  the  pulsed  radar  signal  is 
dehned  as 

N-l 

x{t)  =  J2fkit-kT),  (3.2) 

fc=0 

where 

fkit)  =  afc(t)e^'^'=^*Vect^^^  .  (3.3) 

In  other  words,  the  radar  waveform  is  comprised  of  N  distinct  pulses  separated  by 
a  uniform  interval  T  and  of  uniform  width  Tp.  Subsequently,  when  the  pulse  is 
modulated  onto  a  carrier  it  can  be  denoted  by  the  real  part  of  the  analytic  signal 

frukit)  =  ak{t)  cos(27r/ct  +  0fc(t)).  (3.4) 

Using  the  radar  model.  Figure  2.2,  the  radar  range  channel  r{t)  is  constructed. 
Assuming  the  transmitted  signal  is  received  and  processed  with  coherent  IQ  channels, 
the  analytic  notation  may  be  retained.  Recall  that,  in  the  IQ  representation,  the  real 
part  describes  the  ontpnt  from  the  in-phase  channel,  whilst  the  qnadratnre  channel 
ontpnt  generates  the  imaginary  part.  Substituting  fk(t)  for  x(t)  into  Eq  (2.39) 
describes  the  radar  return  of  the  k^^  pulse, 

r{t)  =  Um{t)-kfk{t) 

=  y  (y )  -  r)  (3.5) 

taking  the  Fonrier  transform  of  both  sides  illustrates  that  the  frequency  response  of 
the  radar  range  channel  is  related  to  the  transmitted  pulse  by 

R{f)  =  Fk{f)Um{f),  (3.6) 
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where 


UM)  = 

=  (3,7) 

which  is  the  spectral  reflectivity  density  of  the  environment,  colored  in  response  to 
an  incident  wave  at  frequency  fc-  Therefore,  the  frequency  response  of  the  radar 
range  channel  is 

R{f)  =  (3.8) 

The  spectrum  of  the  transmitted  pulse  has  a  Altering  effect  on  the  spectral  reflectivity 
density.  High  compression  ratios  from  the  matched  output  predicates  a  requirement 
for  signals  having  large  time-bandwidth  products.  Eq  (3.8)  illustrates  that  large 
bandwidth  signals  also  pass  more  energy  from  the  environment.  The  search  for 
optimal  signals  in  the  Fourier  domain  then  provides  a  more  immediate  relationship 
between  the  signal,  the  radar  return,  and  the  environment. 

3.4-2  Problem  Definition.  Due  to  the  temporal  and  spatial  dynamics  of  the 
spectral  reflectivity  density  function  and  sensor  motion,  each  pulse  return  contains 
different  information  about  the  environment.  The  goal  is  to  extract  this  informa¬ 
tion  from  a  radar  range  channel  with  the  highest  fidelity,  or  minimal  degradation. 
Section  2. 3. 3.1  covers  the  nature  of  range  and  velocity  ambiguities  and  Section  2.4 
discusses  the  benefits  provided  by  NLS.  By  removing  self-clutter  and  ambiguities, 
radar  system  performance  is  improved  and  although  there  are  many  approaches  to 
ambiguity  reduction  the  focus  of  this  thesis  centers  on  the  use  of  a  NLS  technique. 
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Brown  looks  in  the  Fourier  domain  to  divine  an  optimal  solution  which  leads 
to  an  improved  implementation  of  the  NFS  system  [5:18].  Firstly,  complex  notation 
is  used  to  describe  the  spectrum  of  each  of  the  N  pulses  in  the  symbol  set.  Eq  (3.9) 
dehnes  the  Fourier  transform  pair  of  the  pulse 

Fkioo)  =  (3.9) 


where  oj  =  2nf  and  the  envelope  Aiui)  and  phase  <Ffc(i^)  functions  of  Fk{oj)  are 
both  real  and  bandlimited,  he.,  Aioj),  <Ffc(cc;)  G  3?,  with  A{(jj)  =  0  for  all  |  uj\  >  klo/2. 

j: 

However,  whilst  the  signal  Fk{uj)  this  is  not  true  of  the  individual  envelope 

and  phase  function  components,  he.,  A{uj)  ak{t)  and  <hfc(i^)  4>k(t).  The 
components  of  the  NFS  kernel,  in  Figure  3.2,  are  now  redefined  using  this  notation. 
The  matched  filtering  operation  for  the  pulse  is  Hk{uj)  =  F^{uj)  =  A^{(jj)e~AAA ^ 
whilst  the  conjugate  reflection  hlter  is  Gk{uj)  =  [F^{u)]~^  =  A~^  {uj)eAAA _  The 
conjugate  reflection  hlter  from  and  its  neighboring  matched  hlter,  from 
have  reciprocal  envelopes  that  are  independent  of  k.  These  hlters  may  be  combined 
to  form  the  equivalent  hlter  Gk{uj)Hk_i{u)  =  This  hlter  ehectively 

refocuses  the  radar  return,  emphasizing  the  next  pulse,  and  as  such  is  denoted  the 
refocusing  filter.  The  resultant  suppression  kernel,  denoted  the  Brown  kernel  T*,,  is 
dehned  as  illustrated  in  Figure  3.3. 


Figure  3.3.  Brown  threshold  NFS  kernel. 


Note,  also,  that  the  hrst  stage  of  matched  hlters  share  the  envelope  A 2  (a;). 
Subsequently,  the  implementation  of  the  NFS  system  may  be  further  simplihed  by 
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moving  the  envelope  component  A^{(jj)  to  the  front  of  the  suppression  bank.  When 
a  NLS  system  is  constructed  with  the  aforementioned  simplihcations  it  may  be  re¬ 
described  as  shown  in  Figure  3.4.  This  implementation  is  similar  in  structure  to 
the  implementation  shown  in  Figure  2.9  but  is  more  compact  in  terms  of  notation 
and  avoids  singularities  in  the  refocusing  hlter  design.  Therefore,  this  is  the  refer- 


yoit) 

Viif) 


VN-lit) 


Figure  3.4.  Brown  NLS  system  [5]. 


ence  model  used  in  this  thesis  for  implementation  of  Brown  symbols.  With  the  NLS 
system  dehned,  the  problem  becomes  that  of  designing  envelope  A{u)  and  N  phase 
functions  yielding  an  optimal  set  of  symbols.  The  following  section  describes 

Brown’s  approach  to  the  problem. 


3.5  Brown’s  Optimally  Mutual  Dispersive  Symbols 

3.5.1  Optimality  Design  Criteria.  To  provide  a  rigid  mathematical  frame¬ 
work,  compression  and  dispersion  measurements  are  calculated  using  the  appropriate 
metrics  described  in  Section  2. 3. 3. 2.  Recall  from  Eq  (2.49)  that  the  rms  time  dura¬ 
tion  of  a  signal  centered  at  f  =  0  and  having  unit  energy  is  given  by 


a2  = 


0\x{t)\^dt. 


(3,10) 
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Because  of  Fourier  transform  duality,  this  equation  may  be  rewritten  in  terms  of  a 
Fourier  pair  for  time  function  f{t).  In  this  case,  using  the  Fourier  properties 


=  jXiuj) 

1 


\X{u;)\^du; 


(3.11) 

(3.12) 


leads  to  the  equivalent  expression  for  the  rms  time  duration  of  a  signal  in  terms  of 
frequency. 


a2  = 


27r 


jXiu) 


doj. 


(3.13) 


By  dehning  the  spectral  waveform  as  an  arbitrary  complex  analytic  signal,  he., 
X{uj)  =  a  useful  formulation  of  the  rms  time  duration  is  developed: 


a"  = 


1 

1 

1 


—  CXD 

poo 


du 


—  OO 

poo 


1  _ 

P‘^{u)duj  + 


duj 


2ti 


fi^{uj)/3^  {uj)duj. 


(3.14) 


This  is  a  useful  way  of  describing  the  rms  time  duration  because  it  partially  isolates 
the  phase  function,  actually  the  phase-rate  function,  from  the  envelope.  The  NFS 
kernel  performs  thresholding  on  the  correlation  functions.  Consequently,  should  the 
pulse  of  the  symbol  set  dehned  in  Eq  (3.9)  arrive  at  T^.,  then  the  correlation 
function  suppressed  by  the  thresholding  operator  is  dehned  as 


Fi{u;)H,{u;)  =  Fi{u)Fl{u)  =  Pikiu) 

=  (3.15) 


j:  j: 

where  Pik{^)  Pikit)  denotes  the  crosscorrelation  function  and  Pkk{^)  Pkkit) 
denotes  the  autocorrelation  function.  Since  aims  to  suppress  only  the  pulse 
response  the  autocorrelation  function,  pkk(t),  should  be  maximally  compressed  (fo- 
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cused)  whilst  all  crosscorrelation  functions,  pik(t),  should  be  maximally  dispersed. 
The  rms  time  duration  is  small  for  compressed  signals  and  large  for  dispersed  sig¬ 
nals  thereby  providing  a  useful  metric  for  comparing  correlation  functions.  From 
Eq  (3.14),  the  correlation  functions  that  arrive  at  have  the  following  rms  time 
duration 


1  /■“ 

27r 

Where  <Ffc(ci;)  and  A{u)  denotes  the  hrst  derivative  with  respect  to  u.  The  elegance 
of  Brown’s  solution  becomes  no  more  apparent  than  is  shown  here.  For  I  =  k  the  rms 
time  duration  of  the  autocorrelation  function  has  no  phase  component,  leaving  only 
the  first  term  of  Eq  3.16.  The  rms  time  duration  of  the  autocorrelation  function  is 
therefore  independent  of  phase,  decoupling  the  design  of  envelope  A{uj),  to  optimize 
autocorrelation  compression,  from  the  phase  functions  $^(0;),  to  optimize  crosscor¬ 
relation  dispersion.  In  other  words,  rms  time  duration  as  a  metric  of  optimization 
provides  a  means  of  designing  mutually  dispersive  symbols. 


—  ^k{^^)]‘^A'^{u)du.  (3.16) 


1  r°°  . 

afk  =  -  j  ^A\u)du;  + 


3. 5. 1.1  Optimal  Envelope  Design.  Brown  pursues  the  choice  of  an 
envelope  that  optimizes  the  compression  of  a  matched  hltered  bandlimited  signal 
[5:7].  The  design  of  the  optimal  envelope,  A(a;),  such  that  is  minimized  requires 
a  solution  to  the  following  constrained  optimization  problem: 


Minimize 
Such  that 


1 

1 


A‘^{u)duj, 
A^{uj)du  = 


1. 


The  function  achieving  this  is  the  cosine  taper  described  in  Eq  (3.17)  [5:8  . 


A{uj)  = 


0 


,  else 


(3.17) 
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The  cosine  taper,  therefore,  provides  the  optimal  performance  in  rms  time  duration 
of  the  autocorrelation  function  as  illustrated: 


a 


2 

kk 


1 

27r2 

TT^ 


A‘^{uj)duj 


doo 


(3.18) 


The  cosine  taper  is  the  optimal  spectral  envelope  ensuring  minimal  rms  time 
duration  for  a  matched  hltered,  finite  bandwidth  signal.  The  next  step,  therefore, 
is  to  design  phase-rate  functions  providing  the  highest  degree  of  dispersion  in  the 
crosscorrelation  response. 


3. 5. 1.2  Optimal  Phase- Rate  Functions.  The  design  of  optimal  phase- 
rate  functions  begins  by  noticing  that  the  second  term  is  of  Eq  (3.16)  can  be  expanded 
into  components  of  the  form 

{x,y)w=  f  x{t)y*{t)w{t)dt.  (3.19) 

J  a 

This  function  was  hrst  introduced  in  Eq  (2.9)  as  the  weighted  inner  product  of  two 
functions.  To  utilize  this  expression.  Brown  dehnes  the  phase-rate  functions  as  points 
in  a  Hilbert  space,  $^(0;)  G  Ft.  Each  of  the  phase-rate  functions  is  then  dehned  as  a 
linear  combination  of  basis  functions,  (p„(a;),  as  shown 

N-l 

<i>fc(cc;)  =  ^  Ckn  Pn{oo) 

n=0 

=  cl  cpiu;),  (3.20) 
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where 


Cfc  =  [CfcO  Cfcl  Cfc2  •  •  •  CkN-l]'^ 

(p{uj)  =  [ipo{uj)  (pi{u)  (P2{uj)  ... 

G  'H'i  k. 

These  basis  functions  are  chosen  as  odd  functions  of  cn  and  are  orthonormal  within  the 
weighted  inner  product  space  incorporating  A‘^{uj).  The  last  of  these  characteristics 
may  be  expressed  mathematically  as 

=  TT-  /  (pk{uj)(pi{u)A‘^{u)duj  =  5ki  (3.21) 

2vr 

By  incorporating  Eqs  (3.20)  and  (3.21),  the  second  term  of  Eq  (3.16),  describing  the 
dispersion  the  crosscorrelation  function,  is  expanded  as  follows 

2  _  2  _ 

^Ik  ^kk  — 

sub:  dik  =  Cl- Ck 


Subsequently,  the  dispersive  term  in  the  rms  time  duration  can  be  reduced  in  this 
manner  to  the  norm  squared  of  dik,  he.,  the  square  of  the  Euclidean  distance  between 
vectors  c;  and  Ck,  which  is  d\{ci,Ck)  as  dehned  at  Eq  (2.3).  Notably,  by  selecting 
a  satisfactory  basis  set  performance  is  no  longer  dependant  on  the  choice  of  phase- 
rate  functions.  Instead,  the  problem  becomes  one  of  selecting  optimal  coefficient 
weighting  vectors,  he.,  N  vectors  in  iV-dimensional  space  having  the  greatest  equiv- 


1 

^  /  [cf  ‘/^(w)  -  cl  ip{u)]^A\u)du 


1*00  rv  J- 

^  ^  ^lk,n  V^n(^)  ^  ^  ^lk,m  (^UJ^duJ 

n=0  m=0 

A^— 1  A^— 1  y^OO 

EE  dlk,n  dlk^m  I  9^m(^)^  (iS)diO 

^—n  ™ — n  J —oo 


1 

lldzfc 


n=0  m=0 
2 


(3.22) 
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alent  Euclidean  distance.  This  problem  has  been  coined  by  Brown  as  “The  Hermit 
Problem”  [5:9]. 

3.5.2  The  Hermit  Problem.  This  problem  is  so  named  in  reference  to  a 
scenario,  wherein  N-hermits  living  on  a  unit  sphere  in  iV-dimensional  space  jostle  for 
a  position  that  provides  the  greatest  solitude.  The  vectors  addressing  each  hermit  in 
the  steady  state  solution  to  this  competitive  relocation,  being  maximally  equidistant, 
describe  a  solution  to  the  problem.  This  solution,  however,  is  not  unique.  Brown 
begins  his  analysis  of  this  problem  in  2-dimensional  space  [5:9],  eventually  providing 
a  generic  analysis  for  the  more  conceptually  challenging  higher  dimensions. 

Brown  provides  that  in  2-dimensions  the  hermits  will  clearly  reside  at  the 
opposite  ends  of  a  line.  In  3-dimensional  space  they  will  reside  within  the  same 
plane,  each  hermit  separated  by  an  arc  of  120°.  Conceptualization  of  a  solution 
is  confounded  when  considering  4-dimensional  space,  however  Brown  presents  the 
following  generic  process  for  obtaining  a  solution  [5:10].  Begin  with  the  length  N 
vector  a  =  [cos^  cos.^  cos.^  ...  cos.^]^,  where  cos^  =  1/^/N.  This  vector  is  of 
unit  length  and  contains  an  equal  contribution  from  each  of  the  orthogonal  axes. 
Brown  visualizes  this  vector  as  the  handle  of  an  umbrella  with  N  spokes.  As  the 
umbrella  opens,  the  spokes  move  outward  through  the  respective  axes  until  they  all 
become  perpendicular  to  the  umbrella  handle,  a.  Figure  3.5  illustrates  a  particular 
solution  to  the  Hermit  Problem  for  N  =  2  and  N  =  3.  Notably,  a  different  solution 
is  achieved  for  each  different  rotation  of  cc.  However,  this  has  no  impact  on  the 
performance  of  symbols  coded  with  Brown’s  approach. 

Brown  provides  the  following  observations.  The  solution  to  the  hermit  prob¬ 
lem  is  [N  —  l)-dimensional;  in  2-dimensions  the  solutions  exist  on  a  line  and  in 
3-dimensions  within  a  plane,  and  so  on.  Also,  the  sum  of  the  resultant  vectors  is 
zero,  =  0.  Brown  continues,  providing  an  expression  for  optimal  performance 

of  these  symbols  using  a  trigonometric  proof.  The  results  in  Chapter  IV  illustrate 
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1 


Figure  3.5.  A  solution  to  the  Hermit  Problem  for  (a) 
N  =  2,  and  (b)  N  =  3. 


an  implementation  utilizing  a  linear  algebraic  approach,  by  digressing  here  from 
Brown’s  approach  the  solntion  is  obtained  using  a  linear  algebraic  methodology. 

The  nmbrella  analysis  leads  into  the  following  proof  by  emphasizing  the  re- 
qnirement  for  orthogonality.  Each  of  the  solntion  vectors,  c^,  must  be  orthogonal  to 
the  nmbrella  handle  and  as  snch  exist  in  an  orthogonal  (A^— l)-dimensional  space. 
Subseqnently,  the  A^-dimensional  space  may  be  decomposed  into  two  orthogonal 
subspaces.  The  subspace  containing  the  solution  vectors  is  denoted  the  hermit  sub¬ 
space,  (C),  and  the  snbspace  containing  the  nmbrella  handle  is  called  the  orthogonal 
subspace,  (A),  where 


C  =  [cqCi  C2...C7V-1],  (3.23) 

A  =  a.  (3.24) 


The  orthogonal  projection  operator  that  projects  onto  (A)  is  dehned  as  [34:47] 


Pa 


A(A^A)  ^A'^ 

Ct  ^  Ct'^ 


(3.25) 
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where 


1  = 


1  1  1 
1  1  1 
111 


1  1  1  ...  1 


(3.26) 


J  NxN. 


(C)  is  an  orthogonal  subspace  to  (A)  therefore  the  projection  matrix  for  (C) 
is  dehned  as 


'c  = 


Iat  —  Pa 


1 

N 


N  -1 

-1 

-1  ., 

-1 

-1 

N  -1 

-1  ., 

-1 

-1 

-1 

N-l  ., 

-1 

-1 

-1 

-1  ., 

..  iV-1 

(3.27) 
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where  is  the  NxN  identity  matrix.  Each  of  the  N  vectors  in  P(7  is  orthogonal 
to  CK,  also  being  evenly  separated  in  (C).  Subsequently,  a  solution  to  the  Hermit 
Problem  exists  as  the  normalized  column  vectors  of  Pc-  Projection  matrices  are 
idempotent^  therefore  by  observation  the  diagonal  elements  of  Pc  describe  the  norm 
of  the  column  vectors.  Normalizing  the  projection  matrix  Pc  by  dividing  by  the 

^An  idempotent  matrix  when  multiplied  with  itself  results  in  the  original  matrix,  that  is 

n  =Pc. 
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square  root  of  this  value  generates  a  solution  to  the  Hermit  Problem  C  defined  as 


N-l 

-1 

-1  . , 

-1 

r.  -  ^ 

-1 

-1 

N-l 

-1 

-1  . , 

iV-1  ., 

-1 

-1 

y/m-N 

-1 

-1 

-1 

..  iV-1 

(3.28) 


J  NxN. 


The  metric  describing  dispersion  for  Brown’s  approach  is  the  Euclidean  distance. 
This  squared  distance,  between  any  vector  pair  C;,  where  I  ^  k,  is  defined  as 


dl{ci,Ck)  =  dfi^dik 

=  (q  -  Cfc)'^(Cz  -  Cfc) 

2N 

N-l 


(3.29) 


Also,  the  angle  that  is  subtended  between  any  pair  of  vectors  is  given  by 


COS^AT  =  cj  Cfc 
1 

1  -iV 


(3.30) 


Eq  (3.22)  illustrates  that  the  rms  time  duration  for  non-matched,  or  crosscorrelated, 
symbols  is  governed  by  the  Euclidean  distance  between  the  Hermits.  Accordingly, 
the  performance  of  Brown’s  approach  is  provided  by  Eq  (3.29).  This  expression, 
however,  is  given  under  the  constraint  that  each  vector  has  unit  length.  Notably, 
a  uniform  increase  in  the  length  of  the  vectors  will  provide  a  direct  increase  in 
the  Euclidean  distance  withont  compromising  mutuality  and  therefore  produce  an 
improvement  in  the  dispersive  properties  of  the  crosscorrelated  symbols. 
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In  the  following  chapter  dispersive  gain  (Gd),  is  introduced  to  uniformly  in¬ 
crease  the  length  of  the  vector  solutions  to  the  Hermit  Problem,  such  that 

Cfc  ^  \/  G^Ck, 


therefore 


dl{ci,Ck)  =  GD{ci-CkY{ci-Ck) 

=  (^)  Go.  (3.31) 

The  achievable  performance  of  symbols  augmented  by  Gd  is  constrained  only  by 
hardware  limitations. 


3.6  Diserete  NLS  Signals 


3. 6. 1  Desirable  Signal  Properties.  Waveforms  utilized  for  NLS  are  required, 
due  to  the  nature  of  the  threshold  suppression  operator,  to  compress  on  autocorre¬ 
lation  and  disperse  on  crosscorrelation.  Typically,  the  search  for  NLS  codes  has 
leveraged  off  scientihc  and  technical  advances  in  the  held  of  communications,  partic¬ 
ularly  spread-spectrum  communications,  resulting  in  a  range  of  literature  on  the  uses 
of  biphase  and  polyphase  codes.  Such  codes  are  typically  comprised  of  a  sequence 
of  chips,  or  discrete  phase  coded  intervals,  that  concatenate  to  describe  the  symbol, 
or  code.  The  general  form  of  a  complex  phase  coded  signal  is  provided  in  Eq  (3.32) 


[20:152]: 


where 


x{t) 


1 


M-l 

Xm  rect 

m=0 


t  +  -  m) 


X 


m 


0  <  t  <  r 
0,  elsewhere. 


(3.32) 


(3.33) 
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There  are  a  number  of  techniques  used  in  phase-coding  a  sequence  to  achieve  com¬ 
pression,  such  as  Barker  codes  [3]  and  Pseudo-Random  Noise  (PRN)-sequences.  Op¬ 
timally,  these  codes  must  be  mutually  dispersive  to  be  of  any  practical  use  in  NTS. 

3.6.2  Biphase  Codes.  Biphase  codes  are  so  named  because  the  phase- 
code  is  limited  to  being  one  of  two  possible  values,  (pm  G  {0,7r}.  Subsequently,  the 
resultant  code  is  a  real  valued  signal  having  a  chip  value  of  =  ±  1-  These  codes, 
therefore,  are  represented  using  an  M-length  sequence  x  =  e-^*^. 

Barker  sequences  generate  two-valued  autocorrelation  functions,  being  either 
1  or  1/M  [35]  and  thus  have  desirable  autocorrelation  compression  properties.  How¬ 
ever,  only  a  limited  number  of  Barker  codes  exist  and  they  do  not  have  desirable 
crosscorrelation  properties.  Maximal-length  PRN-sequences,  known  as  m-sequences, 
are  usually  generated  via  a  Linear  Feedback  Shift  Register  (LFSR)  [30].  For  m- 
sequences,  the  LFSR  is  conhgured  to  generate  maximal-length  sequences  that  ex¬ 
hibit  two-valued  autocorrelation  functions,  as  per  the  Barker  sequences.  A  Gold 
code  family  is  a  set  of  preferred  pairs  of  m-sequences,  having  crosscorrelation  func¬ 
tions  with  three-values.  Due  to  the  mutually  dispersive  nature  of  Gold  codes,  they 
have  been  employed  in  DS-GDMA  and  spread-spectrum  communications  applica¬ 
tions.  Recently,  these  codes  were  demonstrated  to  provide  an  effective  symbol  set 
for  NLS  [2].  Gold  codes  are  discrete  codes  that  are  well  suited  to  NLS  and  are 
evaluated  in  Section  4.4. 

3. 7  Numerical  Estimation  of  NLS  Symbols  Performance 

3. 7. 1  Performance  Metrics.  The  performance  of  the  symbols,  analyzed  in 
the  following  chapter,  is  measured  in  terms  of  the  rms  time  duration  of  the  correlation 
function,  dehned  in  Eq  (3.34).  The  correlation  function  was  optimized,  in  terms  of 
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rms  time  duration,  for  Brown  symbols  in  Eq  (3.15). 


Pikit) 


/  Ms)fk{s  +  t)ds 

J  — OO 

Fi{u:)F:{u) 


(3.34) 

(3.35) 


The  calculation  of  the  periodic  correlation  function  for  discrete  codes  can  be 
represented  as  shown  in  Eq  (3.36)  [33:594]. 

M-l 

6'n(x,y)  =  A  ^  V  n  =  1-M, -M, . . . ,  M-l,  (3.36) 

m=0 


where 


A 


reel?  ]  dt. 


(3.37) 


"  -  2  '  ^ 

There  are  three  metrics  that  are  commonly  used  to  compare  the  correlation  fnnetions 
for  discrete  codes  [2:3-2].  These  metrics  are:  the  Peak  Autocorrelation  Sidelobe  Level 
(PSL);  Integrated  Autocorrelation  Sidelobe  Level  (ISL);  and  Peak  Crosscorrelation 
Level  (PCCL),  dehned  in  Eqs  (3.38),  (3.39),  and  (3.40)  respectively  [26:537]. 


PSL 


ISL 

PCCL 


10  log 


10  log 

10  log 


I^0(X,X)|2 

■E:::=iMi^™(x,x)n 

m^O 

I^0(X,X)|2 

max  |0„(x,y)|2 

max{|0o(x,x)|2,  |0o(y,y)P} 


(3.38) 


(3.39) 

(3.40) 


The  PSL  metric  provides  a  measure  of  the  disparity  in  the  peak  power  of  the  auto¬ 
correlation  fnnetion  compared  with  the  power  of  the  most  prominent  sidelobe.  ISL 
measures  the  total  power  in  the  sidelobes  compared  to  the  mainlobe.  PCCL  mea¬ 
sures  the  disparity  in  the  peak  power  of  the  antocorrelation  function  compared  with 
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the  power  of  the  most  prominent  sidelobe  in  the  crosscorrelation  function.  There¬ 
fore,  small  sidelobes  will  produce  a  small  PSL,  high  compression  ratios  will  produce  a 
small  ISL  and  high  crosscorrelation  dispersion  will  result  in  a  small  PCCL.  In  short, 
the  smaller  these  values  the  more  suited  the  codes  will  be  to  NLS  implementation. 

Unfortunately,  these  metrics  are  not  dehned  for  continuous  signals.  Therefore, 
the  following  adjustment  is  made  to  accommodate  the  calculation  of  these  param¬ 
eters  for  a  continuous  signal.  Firstly,  the  times  t  at  which  the  extremum  for  the 
autocorrelation  function  occur  must  be  determined,  denoted 

f  =  {t  :  =  0}  .  (3.41) 

The  set,  t,  is  then  used  to  generate  a  discrete  sampling  of  the  autocorrelation  func¬ 
tion,  at  the  extremum,  which  in  turn  leads  to  the  construction  of  the  discrete  auto¬ 
correlation  function 


dn{fkifk)—Pkk{in)i  (3.42) 

where  is  the  ordered  element  of  the  set  t.  This  ordered  set  is  then  used  for 
generating  the  discrete  crosscorrelation  sequences  6{fk,fi).  With  discrete  sequences 
dehned  in  this  way  the  metric  set  is  expanded  to  incorporate  continuous  waveforms. 

3.7.2  Process  Overview.  Chapter  IV  provides  simulation  results  for  vari¬ 
ous  forms  of  Brown  symbols.  The  design  process  divined  from  Section  3.5  may  be 
summarized  as  follows,  for  k  =  0,  1,  . . . ,  iV  — 1: 

1.  Choose  any  set  of  N  unit- vectors  in  iV-dimensional  space  that  solves  The  Her¬ 
mit  Problem.  C  =  [cq  Ci  ...  Cat-i]; 

2.  Choose  any  set  of  N  odd  functions.  (p{uj)  =  [po{oj)  piipj)  . . .  ,  such 

that  ^  jY^ipkY)ipi{u)A^{u)du  =  5ku 
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3.  Describe  the  resultant  N  phase-rate  functions  <hfc(ci;)  =  ^/GdcI 

4.  Describe  the  resultant  N  phase  functions  =  V^dc]^  (p{uj)du-, 

5.  The  resultant  N  mutually  dispersive  codes  expressed  in  the  frequency  domain 

are  then  Fk{ui)  =  and 

6.  The  resultant  N  mutually  dispersive  codes  expressed  in  the  time  domain  are 

Ut)  = 

This  process  is  demonstrated  in  the  following  chapter. 

3.8  Assumptions 

To  meet  the  research  goals  and  stay  within  time  and  resonrce  constraints,  some 
assnmptions  are  necessary.  As  explained  in  Section  3.3,  this  research  seeks  to  simplify 
the  notation  of  Brown’s  approach,  provide  the  theoretical  extension  for  digitization 
of  Brown  symbols,  validate  the  theoretical  expressions  of  rms  time  duration,  and 
analyze  the  resnltant  correlation  waveforms.  As  such  the  following  assumptions 
have  been  made: 

1.  The  Brown  symbols  being  modelled  are  individual  radar  returns  from  station¬ 
ary  single  point  targets; 

2.  The  environment  and  radar  system  do  not  distort  the  phase  information  of  the 
Brown  symbols; 

3.  There  is  no  system  or  environmental  noise;  and 

4.  The  hardware  specihcations  are  snch  that  the  sample  frequency  and  bandwidth 
chosen  for  nnmerical  simulation  are  attainable. 

3. 9  Summary 

This  chapter  described  Brown’s  approach,  which  is  optimal  in  rms  time  dnra- 
tion,  for  generating  a  family  of  mntually  dispersive  symbols  of  any  size.  Section  3.2 
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provided  a  brief  overview  of  the  NLS  concept.  Section  3.3  discussed  the  scope  of 
the  research,  relating  the  approach  and  associated  synthesis  constraints.  Section  3.4 
related  the  radar  model  to  NLS  and  analytically  defined,  as  in  [5],  the  problem  of 
designing  a  symbol  set  suited  to  NLS.  Section  3.5  derived  Brown’s  theorem  divining 
the  optimal  envelope,  phase-rate  functions  and  solutions  to  the  Hermit  Problem,  and 
also  introduced  dispersive  gain  and  provided  expressions  for  optimal  performance. 
Section  3.6  dehned  notation  for  biphase  codes  relating  the  suitability  of  these  codes 
to  NLS.  Finally,  Section  3.7  discussed  the  metrics  used  to  characterize  and  contrast 
Brown  symbols,  providing  an  overview  of  the  numerical  approximation  approach 
undertaken. 
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IV.  Analysis  and  Results 


4-1  Introduction 

In  this  chapter,  Brown  symbols  are  devised  analytically,  using  the  process 
outlined  in  Section  3.7.2,  and  then  constructed  and  simulated  in  Matlab®.  Section  4.2 
considers  the  design  of  each  of  the  Brown  symbol  components,  being  the  Hermit 
Problem  solution,  phase-rate  basis  function  design,  phase  function  calculation,  and 
symbol  synthesis.  Section  4.3  demonstrates  the  invariance  of  Brown  symbols  to  the 
choice  of  phase-rate  basis  functions  and  envelope,  and  considers  physical  realizability. 
Finally,  Section  4.4  contrasts  the  performance  of  Brown  symbols  with  that  of  extant 
biphase  symbols. 


4-2  Brown  Symbol  Design 

Section  3. 5. 1.1  devised  the  optimal  envelope  for  Brown  symbols.  This  section 
describes  the  process  of  generating  Brown  symbols  by  solving  the  Hermit  Problem 
and  designing  phase  functions. 


4-2.1  Hermit  Problem  Solution.  Section  3.5.2  outlined  a  methodology 
for  solving  the  Hermit  Problem  by  using  a  linear  algebraic  technique.  Given  the 
starting  vector  a  described  in  Section  3.5.2,  this  approach  provides  a  simple  means 
of  generating  the  solution  matrix,  C,  as  detailed  in  Eq  (4.1). 


N1n-1 


(4.1) 
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where  In  is  an  N  x  N  identity  matrix,  and 


1  1  1  ...  1 

1  1  1  ...  1 


(4.2) 


1  1  1 


1 


NxN. 


Alternative  solutions  to  the  Hermit  Problem  can  be  attained  by  reworking  the 
derivation  using  a  different  vector  for  the  umbrella  handel,  a,  or  by  operating  on 
the  solution  provided  here  with  an  iV-dimensional  rotation  matrix.  Once  a  solution 
has  been  calculated,  the  next  step  of  generating  Brown  symbols  is  to  determine  the 
phase  functions. 


4-2.2  Phase  Function  Design.  Section  3.5. 1.2  provided  a  methodology  for 
the  design  of  optimal  phase-rate  functions  and  indicated  that  N  basis  functions  are 
required.  These  functions  must  be  odd  functions  of  u  and  orthonormal  within  the 
weighted  inner  product  space  incorporating  A^(a;).  A  generic  four-step  process  for 
generating  these  basis  functions  follows: 

1.  Choose  any  set  of  N  odd  and  linearly  independent  functions. 

'd{u)  = 

2.  Using  any  orthonormalization  procedure  with  the  weighted  inner-product,  gen¬ 
erate  a  set  of  orthonormal  basis  functions  ■  ■  ■  Pn-i{!.jj)Y  - 

3.  Describe  the  resultant  N  phase-rate  functions  where 

Gd  is  the  dispersive  gain  factor  and  is  the  solution  to  the  Hermit  Prob¬ 
lem. 

4.  Calculate  the  resultant  N  phase  functions  =  VGdY:  fuj  p{^)duj. 
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This  process  can  be  very  involved  but  permits  the  widest  selection  of  basis  functions 
for  implementation  in  Brown  symbols.  Once  the  phase  functions  are  calculated  the 
Brown  symbols  can  be  compiled. 


4-2.3  Symbol  Synthesis.  The  phase  functions  are  the  key  element  in  de¬ 
signing  strong  mutually  dispersive  symbols.  Once  the  N  phase  functions  have  been 
calculated  the  Brown  symbols  maybe  compiled.  The  resultant  iV-symbol  family  is 
represented,  for  A;  =  0,  1,  . . . ,  iV,  as 


Fk{uj)  = 


fk{t)  = 


(4.3) 

(4.4) 


The  phase  functions  are  sampled  and  reconstructed  as  vectors  in  Matlab®.  Care  is 
taken  to  ensure  that  the  frequency  sample  interval  is  chosen  to  avoid  time-aliasing. 
The  element  of  the  sampled  phase  function  is  described  as  ^kn  = 
the  tilde  notation  is  used  to  denote  the  sampled  approximation  to  the  continuous 
frequency  or  time  function  .  Due  to  sampling  in  the  frequency  domain  the  corre¬ 
sponding  time  function  will  be  periodic  with  period  Tq  =  ‘I'n therefore  loq  is 
chosen  to  avoid  aliasing  such  that  loq  <  /Tq.  This  is  simply  the  converse  to  the 
time  sampled  analysis  provided  in  Section  2.2.3.  Each  of  the  phase  functions  and 
the  envelope  A„  =  A{nujQ)  are  sampled  in  this  manner  to  synthesize  the  spectral 
Brown  symbols.  The  element  of  the  Brown  symbol,  in  frequency,  is  therefore 
dehned  as 


Ykn  =  Ke^^^-.  (4.5) 

As  mentioned  earlier,  these  symbols  are  constructed  in  Matlab®.  Matlab®  is  a  dis¬ 
crete  simulation  application  that  constrains  further  analysis  to  a  numerical  approx¬ 
imation.  The  corresponding  time  domain  signal  for  each  of  the  Brown  symbols  is 
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calculated  using  the  rectangular-rule  approximation  to  the  Fourier  transform,  which 
is  the  DFT  as  shown  in  Eq  (2.33).  The  resultant  vector  has  an  equivalent  number 
of  samples  (iV^)  and  a  time  sample  interval  to  =  ^tt/NsUq.  The  element  of  the 
^th  symbol,  in  time,  is  dehned  as 


f 


kn 


1 

N. 


Ns-l 


—  Fki 


1 


1=0 

Ns-l 


A;  Q^\-^kl+‘iTrln/Ns 


1=0 


(4.6) 


Consequently,  the  results  provided  in  this  chapter  are  parameterized  by  the 
following:  choice  of  basis  number  of  symbols  (iV);  bandwidth  (hlo);  and 

dispersive  gain  (Gd).  The  last  of  these,  dispersive  gain,  was  introduced  in  Section 
3.5.2  as  a  factor  that  uniformly  increases  the  length  of  the  vectors  that  solve  the 
Hermit  Problem,  thereby  controlling  the  dispersion  induced  in  the  crosscorrelation 
functions. 


4-3  Brown  Symbol  Sensitivity 

Brown’s  theorem  indicates  that  the  resultant  symbols  exhibit  properties  that 
are  insensitive  to  changes  in  basis.  However,  this  is  only  in  consideration  of  the  rms 
time  duration  of  the  correlation  functions.  This  section  considers  a  number  of  differ¬ 
ent  bases  illustrating  that  although  the  theoretical  results  are  achieved  independent 
of  basis,  the  form  of  the  correlation  functions  vary. 

Unless  otherwise  stipulated,  the  results  provided  in  this  section  are  generated 
using  a  bandwidth  of  hlo  =  Under  this  condition  the  theoretical  rms  time 
duration  for  the  autocorrelation  function  is  =  0.0156,  and  for  the 

crosscorrelation  functions  is  =  Tr^/hlo  +  {‘^GdN)/{N  —  1),  which  for  a  symbol 
family  of  iV  =  3,  and  Gn  =  0  dB,  is  3.0156. 
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4-3.1  Sinusoidal  Basis.  Sinusoids  are  an  obvious  choice  as  basis  functions 
due  to  their  orthogonality  properties.  Appendix  B.l  derives  the  generation  of  an 
A^-dimensional  set  of  sinusoidal  basis  functions  that  retain  orthonormality  within 
the  weighted  inner  product  space.  The  general  form  for  each  element  in  this  set  of 
basis  functions  is  described  in  Eq  (4.7). 

(fniuj)  =  y2sin(ano;),  (4.7) 


therefore 


*hfc(cc;) 


N-l 


E\/2G£)Ckn  ,  N 

- cos(ano;), 


n=0 


(4.8) 


where  a„  is  an  arbitrary  positive  integer  and  c^^  is  the  element  of  the  solution 
to  the  Hermit  Problem.  As  indicated  in  Appendix  B.l,  for  orthogonality  to  be 
exactly  met,  this  basis  requires  the  bandwidth  to  be  an  even  integer  multiple  of  tt, 
he.,  Hq  =  2m7r,  m  E  Z.  The  synthesized  Brown  symbols  are  then  of  the  form 


Fk{uj) 


fkit) 


COs(an<.^) 


3 

e 


N-l 

n=0 


V^GpCfc,,  cos(a„a4) 

Cln 


dui. 


(4.9) 

(4.10) 


A  closed  form  solution  to  Eq  (4.10)  is  unlikely  to  exist,  however,  a  numerical  ap¬ 
proximation  to  this  function  provides  a  means  of  analysis.  This  approach  does 
not  detract  from  the  end  result  of  physical  implementation  as  the  sampled  form  of 
fk(t)  tofk{nto)  =  fkn  is  the  code  that  a  digital  radar  system  would  utilize,  with 
parameters  such  as  sample  interval  governed  by  hardware  specihcs. 

These  symbols  are  synthesized,  as  detailed  in  Section  4.2.3,  for  a  code  family 
of  A^  =  3  symbols.  A  dB  plot  of  the  autocorrelation  function,  which  is  the  same  for 
each  symbol,  is  displayed  in  Eigure  4.1.  The  x-axis  describes  the  time  lag  and  the 


4-5 


y-axis  describes  the  amplitude,  in  decibels.  The  metrics  described  in  Table  4.1  have 


Figure  4.1.  Autocorrelation  function  for  a  sinusoidal  basis  with 
AT  =  3,  Gd  =  0  dB,  and  a  =  [13  18  23]^. 

been  calculated  from  these  crosscorrelation  functions. 

Table  4.1.  Autocorrelation  metrics  for  a  sinusoidal  basis  with  N  =  3, 
Gd  =  0  dB,  and  a  =  [13  18  23]^. 


Plk{t) 

PSL 

IdBI 

ISL 

ldB| 

Poo{i) 

-23.000 

-18.931 

0.016 

Piiii) 

-23.000 

-18.931 

0.016 

P22{t) 

-23.000 

-18.931 

0.016 

A  dB  plot  of  the  crosscorrelation  functions  is  displayed  in  Figure  4.2.  The 
x-axis  describes  the  time  lag  and  the  y-axis  describes  the  amplitude,  in  decibels. 
The  metrics  described  in  Table  4.2  have  been  calculated  from  these  crosscorrelation 
functions.  Although  the  theoretical  results  for  and  ah  are  achieved  the  cross¬ 
correlation  plots  do  not  show  strong  dispersive  properties  resulting  in  undesirably 
small  values  of  PCCL  and  dominant  sidelobes.  By  increasing  the  dispersive  gain 
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Figure  4.2.  Crosscorrelation  functions  for  a  sinusoidal  basis  with 
N  =  3,  Gd  =  0  dB,  and  a  =  [13  18  23]'^. 


Table  4.2.  Crosscorrelation  metrics  for  a  sinusoidal  basis  with  N  =  3, 
Gd  =  0  dB,  and  a  =  [13  18  23]^. 


Plk{t) 

PCCL 

|dB] 

Poi{t) 

-0.059 

3.016 

P02{t) 

-0.051 

3.016 

Puit) 

-0.032 

3.016 

4-7 


the  amount  of  spread  is  controlled.  Figures  4.3  and  4.4  illustrate  the  effect  on  the 
crosscorrelation  functions  of  increasing  Go  to  20  and  40  dB  respectively.  The 
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Figure  4.3.  Crosscorrelation  functions  for  a  sinusoidal  basis  with 
N  =  3,  Gd  =  20  dB,  and  a  =  [13  18  23]^. 


Table  4.3.  Crosscorrelation  metrics  for  a  sinusoidal  basis  with  N  =  3, 
Gd  =  20  dB,  and  a  =  [13  18  23]^. 


Plkit) 

PCCL 

[dB] 

Poi{t) 

-6.537 

300.016 

P02{t) 

-5.670 

300.016 

Pl2{t) 

-3.420 

300.016 

metrics  described  in  Tables  4.3  and  4.4  have  been  calculated  from  these  correlation 
functions  shown  in  Figures  4.3  and  4.4  respectively. 

Increasing  the  dispersive  gain.  Go,  has  the  effect  of  improving  the  dispersive 
properties  of  the  crosscorrelation  functions  without  affecting  the  autocorrelation 
functions.  To  improve  the  compression  of  the  autocorrelation  function  the  band¬ 
width  must  be  increased  as  illustrated  (for  a  sinusoidal  basis)  in  Figure  4.5.  This 
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Figure  4.4.  Crosscorrelation  functions  for  a  sinusoidal  basis  with 
AT  =  3,  Gd  =  40  dB,  and  a  =  [13  18  23]^. 


Table  4.4.  Crosscorrelation  metrics  for  a  sinusoidal  basis  with  iV  =  3, 
Gd  =  40  dB,  and  a  =  [13  18  23]^. 


Plkit) 

PCCL 

[dB] 

(^ki 

Pm{t) 

-19.985 

30000.016 

P02{t) 

-14.925 

30000.016 

Pl2{t) 

-18.789 

30000.016 
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basis  achieves  the  theoretical  Brown  autocorrelation  performance,  as  a  function  of 
bandwidth. 


0.25 


0.2 

0.15 

3 

0.1 

0.05 


Figure  4.5.  for  a  sinusoidal  basis  with  Gd  =  0  dB  and 

a  =  [13  18  23]^. 

The  performance  of  the  crosscorrelation  in  terms  of  rms  time  duration  is  con¬ 
sidered  in  relation  to  the  number  of  symbols  (iV).  Figure  4.6  indicates  that  the 
theoretical  Brown  crosscorrelation  result,  for  symbol  families  with  N  =  2,  3,  . . . ,  7, 
is  also  achieved  with  the  sinusoidal  basis.  Although  the  crosscorrelation  dispersion 
decreases  with  an  increase  in  the  number  of  symbols,  increasing  the  dispersive  gain 
offsets  this  effect. 

This  section  demonstrated  that  the  sinusoidal  basis  is  a  valid  basis  for  any  size 
symbol  family,  achieving  the  theoretical  Brown  results  in  each  case.  This,  however, 
says  nothing  of  the  suitability  of  these  symbols  for  physical  implementation  or  their 
comparative  performance.  These  aspects  are  now  considered. 

4 .3. 1.1  Realizable  Brown  Symbols.  Brown  symbols  are  bandlimited 
and  therefore  unlimited  in  time.  However,  a  physical  implementation  of  these  sym- 
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Figure  4.6.  —  ah.  for  a  sinusoidal  basis  with  Gd  =  0  dB  and 

a  =  [13  18  23]^. 

bols  requires  them  to  be  time-limited.  Consequently,  an  analysis  of  windowed  Brown 
symbols  is  provided.  Following  this  a  comparison  is  made  between  the  performance 
of  a  pair  of  sinusoidal  Brown  symbols  and  the  conjugate  LFM  symbols  used  in  the 
initial  demonstration  by  Palermo  [28]. 

The  requirement  to  transmit  time-limited  symbols  leads  to  an  analysis  of  win¬ 
dowed  Brown  symbols.  In  this  analysis  a  rectangular  window,  w{t),  of  unit  height  is 
employed  to  create  a  time-limited  Brown  symbol, 

n;(t)  =  rect(^^^  ,  (4.11) 

where  is  the  width  of  the  window.  The  window  will  zero  the  tails  of  the  Brown 
symbols  for  fk(t),  \t\  >  Ty^/2.  Figure  4.7  illustrates,  for  Brown  symbols  generated 
with  various  dispersive  gains,  that  the  signal  energy  becomes  more  concentrated  as 
the  dispersive  gain  is  increased.  The  x-axis  describes  the  dispersive  gain,  in  decibels. 
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and  the  y-axis  describes  the  relative  half-window  width  Tu,/2  as  a  fraction  of  the  rms 
time  duration  of  the  Brown  symbol  cifc. 


Figure  4.7.  Relative  window  duration  T^/2ak  required  to  capture 
99.9%  of  signal  energy  as  a  function  of  dispersive  gain 

Gd- 

Accordingly,  Brown  symbols  generated  using  a  large  dispersive  gain  are  more 
compact  and  are  therefore  better  suited  to  windowing.  Figure  4.8  displays  the  ef¬ 
fect  that  shortening  the  window  width  has  on  the  properties  of  the  crosscorrelation 
function.  Notably,  the  effect  of  windowing  Brown  symbols  becomes  insignihcant  for 
Tw/2  K,  2.5  (Tfc.  Therefore  the  performance  of  Brown  symbols  is  not  signihcantly 
degraded  due  to  windowing,  for  large  dispersive  gain,  suggesting  that  these  symbols 
are  easily  modihed  for  physical  implementation. 

Originally  Palermo  utilized  conjugate  LFM  pulses  to  demonstrate  the  effec¬ 
tiveness  of  NFS.  Unfortunately  the  number  of  conjugate  LFM  is  limited  to  iV  =  2. 
Having  demonstrated  that  Brown  symbols  may  be  windowed  for  physical  imple¬ 
mentation  a  comparison  is  now  provided  between  the  characteristics  of  sinusoidal 
Brown  symbols  and  conjugate  LFM.  In  this  comparison  symbols  are  generated  hav- 
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Figure  4.8.  Correlation  function  statistics  for  windowed  Brown  sym¬ 
bols  with  sinusoidal  basis  for  N  =  2,  Gu  =  40  dB,  and 
a  =  [13  18]^.  (a)  PSL,  (b)  auk,  (c)  PCCL,  and  (d)  an- 


ing  equivalent  rms  time  duration  and  energy.  The  results,  in  Figure  4.9,  illustrate 
the  comparative  performance.  Brown  symbols,  being  optimized  for  rms  time  dura¬ 
tion  of  the  correlation  function,  show  equivalent  or  better  performance  in  all  areas 
except  the  peak  crosscorrelation  level  (PCCL).  Figure  4.9  (c)  shows  that  the  Brown 
symbols  exhibit  approximately  10  dB  worse  PCCL  than  conjugate  LFM.  However, 
Brown  symbols  exhibit  consistent  performance  when  compared  with  conjugate  LFM, 
and  unlike  conjugate  LFM,  Brown  symbols  are  not  limited  in  symbol  family  size. 


4-3.2  Polynomial  Basis.  Polynomials  have  been  shown  in  various  applica¬ 
tions  to  exhibit  strong  orthogonality  properties,  as  illustrated  by  the  set  of  Chebyshev 
polynomials,  Hermite  functions,  Legendre  polynomials  and  so  on  [36:136-160]  [16:25- 
36].  Unfortunately,  these  specihc  sets  are  orthogonal  only  within  a  unity  weighted 
inner  product  space  requiring  further  analysis  for  adaptation  to  optimal  Brown  sym¬ 
bols.  The  design  of  polynomial  basis  for  use  in  Brown  codes  is  explored  in  this 
section. 
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Figure  4.9.  Correlation  function  statistics  for  conjugate  LFM  and 
Brown  symbols  with  sinusoidal  basis  for  N  =  2  and 
a  =  [13  18]^.  (a)  PSL,  (b)  akk,  (c)  PCCL,  and  (d) 

O'kl- 


Appendix  B.2  shows  that  an  iV- dimensional  set  of  polynomial  basis  functions, 
which  retain  orthonormality  within  the  weighted  inner  product  space,  can  be  gener¬ 
ated.  The  general  form  for  each  element  in  this  set  of  basis  functions,  is  described 
in  Eq  (4.12).  Each  basis  element  of  order  2p  —  1,  is  a  weighted  set  of  odd 

polynomials: 

p-i 

V9n(cn)  =  (4.12) 

m=0 

where  hnm  is  the  element  in  the  row  of  the  column  of  the  matrix  B  dehned 
in  Eq  (B.12).  Therefore 


N-l 


p-1 


*hfc(ci..’)  ^  ^  's/ G £)Ckn  ^  ^  bf; 


00 


2m+2 


n=0 


m=0 


‘2m  +  2 


=  f  ip{oj)doj, 

J  u 


(4.13) 
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where  is  the  solution  to  the  Hermit  Problem.  The  synthesized  Brown  symbols 
are  therefore 


Fki,^) 

fkit) 


recti  — 
\Iq 


—  cos5  (  —  I 


Hfi 


duj. 


(4.14) 

(4.15) 


As  was  the  case  with  Eq  (4.10),  a  closed  form  solution  to  Eq  (4.15)  is  unlikely  to 
exist,  however,  a  numerical  approximation  to  this  function  again  provides  a  means 
of  analysis. 

These  symbols  are  synthesized,  as  detailed  in  Section  4.2.3,  for  a  code  family 
of  TV  =  3  symbols.  A  dB  plot  of  the  antocorrelation  fnnction,  which  is  the  same  for 
each  symbol,  is  displayed  in  Figure  4.10.  The  x-axis  describes  the  time  lag  and  the 
y-axis  describes  the  amplitnde,  in  decibels.  The  metrics  described  in  Table  4.5  have 


Fignre  4.10.  Antocorrelation  fnnction  for  a  polynomial  basis  with 
TV  =  3  and  Gd  =  ^  dB. 

been  calculated  from  these  crosscorrelation  fnnctions. 
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Table  4.5.  Autocorrelation  metrics  for  a  polynomial  basis  with  N  =  3 
and  Gd  =  0  dB. 


Plk{t) 

PSL 

ldB| 

ISL 

ldB| 

Poo(t) 

-23.048 

-18.975 

0.016 

Pll(t) 

-23.048 

-18.975 

0.016 

P22{t) 

-23.048 

-18.975 

0.016 

A  dB  plot  of  the  crosscorrelation  functions  is  displayed  in  Figure  4.11.  The 
x-axis  describes  the  time  lag  and  the  y-axis  describes  the  amplitude,  in  decibels.  The 


Figure  4.11.  Crosscorrelation  functions  for  a  polynomial  basis  with 
N  =  3  and  Gd  =  0  dB. 

metrics  described  in  Table  4.6  have  been  calculated  from  these  correlation  functions. 


Again,  although  the  polynomial  basis  has  achieved  the  theoretical  results  for 
all,  these  crosscorrelation  plots  also  fail  to  show  strong  dispersive  properties. 

By  increasing  the  dispersive  gain  the  amount  of  spread  is  controlled.  Figures  4.12 
and  4.13  illustrate  the  effect  on  the  crosscorrelation  functions  of  increasing  Gd  to 
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Table  4.6.  Crosscorrelation  metrics  for  a  polynomial  basis  with  N  =  3 
and  Gd  =  0  dB. 


Plk{t) 

PCCL 

|dB] 

PoiW 

-6.482 

3.016 

P02{t) 

-7.690 

3.016 

Pl2{t) 

-7.382 

3.016 

20  and  40  dB  respectively.  The  metrics  described  in  Tables  4.7  and  4.8  have 
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Figure  4.12.  Crosscorrelation  functions  for  a  polynomial  basis  with 
N  =  3  and  =  20  dB. 

been  calculated  from  the  crosscorrelation  functions  shown  in  Figures  4.12  and  4.13 
respectively.  This  section  demonstrated  that  the  polynomial  basis  is  a  valid  basis 
for  any  size  symbol  family,  achieving  the  theoretical  results  in  each  case. 

4-3.3  Chebyshev  Basis.  Appendix  B.3  shows  that  an  A^-dimensional  set 
of  Chebyshev  basis  functions  can  be  generated  when  a  non-ideal  uniform  spectral 
envelope  is  employed.  Under  this  condition  the  element  in  the  set  of  Chebyshev 


4-17 


0 


-300  -200  -100  0  100  200  300 


-500  -400  -300  -200  -100  0  100  200  300  400  500 


Figure  4.13.  Crosscorrelation  functions  for  a  polynomial  basis  with 
N  =  3  and  Gd  =  40  dB. 


Table  4.7.  Crosscorrelation  metrics  for  a  polynomial  basis  with  N  =  3 
and  Gd  =  20  dB. 


Plk{t) 

PCCL 

[dB] 

(^ki 

Poi{t) 

-15.407 

300.016 

P02{t) 

-15.187 

300.016 

Pl2{t) 

-11.159 

300.016 

Table  4.8.  Crosscorrelation  metrics  for  a  polynomial  basis  with  N  =  3 
and  Gd  =  40  dB. 


Pikit) 

PCCL 

[dB] 

Pm{t) 

-22.189 

30000.016 

P02{t) 

-22.264 

30000.016 

Pl2it) 

-15.978 

30000.016 
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basis  functions  is  described  in  Eq  (4.16) 


‘-Pniuj) 


A 


(4.16) 


where  Tfe(ci;)  is  a  order  Chebyshev  polynomial  defined  in  Eq  (B.14).  For  or¬ 
thonormality  the  envelope  is  uniform  over  flo  and  is  dehned  as 


and  the  phase  function  of  the  symbol  is  dehned  as 


(4.17) 


AT-l 


E 

n=0 


(4.18) 


where  Ckn  is  the  element  of  the  k^^  solution  to  the  Hermit  Problem.  The  synthe¬ 
sized  Brown  symbols  are  then 


Fk{uj) 

fk{t) 


(4.19) 

(4.20) 


Again,  a  closed  form  solution  to  Eq  (4.20)  is  unlikely  to  exist,  however,  a  numerical 
approximation  to  this  function  provides  a  means  of  analysis. 

These  symbols  are  synthesized,  as  detailed  in  Section  4.2.3,  for  a  code  family 
of  TV  =  3  symbols.  A  dB  plot  of  the  autocorrelation  function,  which  is  the  same  for 
each  symbol,  is  displayed  in  Figure  4.14.  The  x-axis  describes  the  time  lag  and  the 
y-axis  describes  the  amplitude,  in  decibels.  The  metrics  described  in  Table  4.9  have 
been  calculated  from  these  crosscorrelation  functions. 
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Figure  4.14. 


Autocorrelation  function  for  a  Chebyshev  basis  with 
N  =  3  and  Gn  =  0  dB. 


Table  4.9.  Autocorrelation  metrics  for  a  Chebyshev  basis  with  N  =  3 
and  Gd  =  0  dB. 


Plkit) 

PSL 

[dB| 

ISL 

|dB] 

^kk 

Pooit) 

-13.267 

-7.076 

316.033 

Pll{t) 

-13.267 

-7.076 

316.033 

P22{t) 

-13.267 

-7.076 

316.033 
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A  dB  plot  of  the  crosscorrelation  functions  is  displayed  in  Figure  4.15.  The 
x-axis  describes  the  time  lag  and  the  y-axis  describes  the  amplitude,  in  decibels. 
The  metrics  described  in  Table  4.10  have  been  calculated  from  these  crosscorrelation 
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Figure  4.15.  Crosscorrelation  functions  for  a  Chebyshev  basis  with 
N  =  3  and  Gd  =  0  dB. 

functions. 

Table  4.10.  Crosscorrelation  metrics  for  a  Chebyshev  basis  with 
N  =  3  and  Gd  =  0  dB. 


plkit) 

PCCL 

[dB] 

Poi(t) 

-8.983 

319.033 

P02{t) 

-8.564 

319.033 

Puit) 

-9.034 

319.033 

Although  the  Chebyshev  basis  provides  worse  autocorrelation  compression,  for 
the  same  dispersive  gain  the  PCCLs  are  better  than  for  the  sinusoidal  and  poly¬ 
nomial  bases.  Notably,  a  the  non-ideal  envelope  provides  worse  compression  but 
better  dispersion,  normally  this  would  lead  to  an  envelope  design  trade-off  between 
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dispersion  and  compression.  However,  the  dispersive  gain  controls  dispersion  enough 
to  accommodate  for  any  trade-off  that  might  otherwise  have  been  necessary.  Fig¬ 
ures  4.16  and  4.17  illustrate  the  effect  on  the  crosscorrelation  functions  of  increasing 
Gd  to  20  and  40  dB  respectively.  The  metrics  described  in  Tables  4.11  and  4.12 


Figure  4.16.  Crosscorrelation  functions  for  a  Chebyshev  basis  with 
N  =  3  and  Gd  =  20  dB. 


Table  4.11.  Crosscorrelation  metrics  for  a  Chebyshev  basis  with 
N  =  3  and  Gd  =  20  dB. 


Plk{t) 

PCCL 

[dB] 

Pm{t) 

-16.066 

616.029 

P02{t) 

-15.617 

616.029 

Pl2{t) 

-16.830 

616.029 

have  been  calculated  from  these  correlation  functions  shown  in  Figures  4.16  and  4.17 
respectively. 

As  expected,  the  change  in  A{uj)  to  a  non-ideal  rectangular  window  results  in 
worse  performance  of  the  autocorrelation  function,  he.,  larger  however  mutual 
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Figure  4.17.  Crosscorrelation  functions  for  a  Chebyshev  basis  with 
N  =  3  and  =  40  dB. 


Table  4.12.  Crosscorrelation  metrics  for  a  Chebyshev  basis  with 
N  =  3  and  Gd  =  40  dB. 


Plkit) 

PCCL 

[dB] 

(^ki 

Poi{t) 

-22.780 

30315.616 

P02{t) 

-22.655 

30315.616 

Pl2{t) 

-24.680 

30315.616 
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dispersion  is  exhibited  due  to  orthonormality  in  the  basis  set.  The  control  on  dis¬ 
persion  offered  by  the  dispersive  gain  accommodates  any  envelope  design  trade-off 
that  might  be  made  of  worse  compression  for  increased  dispersion.  This  section 
analyzes  the  Chebyshev  basis,  for  any  size  symbol  family,  showing  that  although 
the  results  are  suboptimal  the  characteristics  of  the  crosscorrelation  functions  are 
mutually  dispersive  and  well  suited  to  NLS. 


4.3.4  Piecewise  Continuous  Basis.  Appendix  B.4  shows  that  a  set  of  N 
orthonormal  piecewise  basis  functions  can  be  generated  using  the  uniform  taper 
applied  previously  for  the  Chebyshev  case.  The  general  form  for  each  element  in  this 
set  of  basis  functions  is  dehned  in  Eq  (B.19)  as 


.  .  A  4^/3Sn 
Pn{uj)  =  - 

23^71  1 


I  I  n  /o^  [ W)  +  0.^12  — 

(a,  +  SJ„/2)rect(^ - - 


+  (-1)”‘(^  +  u',ii)rect(  ^  ^ 

m=l  ^  ' 

<  t  n  /o^  ~  ^0/2  + 

+  (a,  -  n„/2)rect(^ - - 


(4.21) 


therefore 


N-l 


V Gd> 


D^kn " 


n=0 
2371  —  1 


A^/3c 


+  r2oi^)rect' 


CO  r2o/2  —  C,n/8 


m= 

,2 


+  2o;^o;)rect(  ^7^" 

m=l 

oj  —  VLq/2  +  VLn/8 


+  —  r2oi^)rect 


ffj4 


(4.22) 


where  Ckn  is  the  element  of  the  solution  to  the  Hermit  Problem.  The  synthe¬ 
sized  Brown  symbols  are  then  of  the  form 


Pk{uj)  = 


e 


(4.23) 
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1 


(4.24) 


fk{t)  = 


TTOO 

Qq 


Again,  a  closed  form  solution  to  Eq  (4.24)  is  unlikely  to  exist,  so  a  numerical  ap¬ 
proximation  to  this  function  is  considered. 

These  symbols  are  synthesized,  as  detailed  in  Section  4.2.3,  for  a  code  family 
of  TV  =  3  symbols.  A  dB  plot  of  the  autocorrelation  function,  which  is  the  same 
for  each  symbol,  is  displayed  in  Figure  4.18.  The  x-axis  describes  the  time  lag  and 
the  y-axis  describes  the  amplitude,  in  decibels.  The  metrics  described  in  Table  4.13 


Figure  4.18.  Autocorrelation  function  for  a  piecewise  basis  with 
N  =  3,  Gd  =  0  dB,  and  a  =  [12  A]'^. 

have  been  calculated  from  these  crosscorrelation  functions. 

A  dB  plot  of  the  crosscorrelation  functions  is  displayed  in  Figure  4.19.  The 
x-axis  describes  the  time  lag  and  the  y-axis  describes  the  amplitude,  in  decibels. 
The  metrics  described  in  Table  4.14  have  been  calculated  from  these  crosscorrelation 
functions.  The  correlation  plots  exhibit  similar  characteristics  to  the  Chebyshev 
basis,  having  less  compression  on  the  autocorrelation  function  and  slightly  better 
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Table  4.13.  Autocorrelation  metrics  for  a  piecewise  basis  with  TV  =  3, 
Gd  =  0  dB,  and  a  =  [1  2  4]^. 


Plk{t) 

PSL 

ldB| 

ISL 

|dB| 

Pooit) 

-13.262 

-7.109 

63.203 

Pll{t) 

-13.262 

-7.109 

63.203 

P22{t) 

-13.262 

-7.109 

63.203 
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Figure  4.19.  Crosscorrelation  functions  for  a  piecewise  basis  with 
iV  =  3,  Gd  =  0  dB,  and  a  =  [12  4]^. 


Table  4.14.  Crosscorrelation  metrics  for  a  piecewise  basis  with  iV  =  3, 
Gd  =  0  dB,  and  a  =  [12  4]'^. 


Pikit) 

PCCL 

|dB] 

Poi{t) 

-8.187 

66.203 

P02{t) 

-9.307 

66.203 

Pl2{t) 

-4.788 

66.203 
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dispersion  on  the  crosscorrelation  functions.  The  ort honor mality  in  the  basis  pro¬ 
vides  the  theoretical  values  for  the  dispersive  term  and  the  form  of  the  correlation 
functions  is  suited  to  NLS,  being  similar  to  that  of  the  Chebyshev  and  polynomial 
bases.  Figures  4.20  and  4.21  illustrate  the  effect  on  the  crosscorrelation  functions  of 
increasing  Gd  to  20  and  40  dB  respectively.  The  metrics  described  in  Tables  4.15 


Figure  4.20.  Crosscorrelation  functions  for  a  piecewise  basis  with 
N  =  3,  Gd  =  20  dB,  and  a  =  [12  A]'^. 


Table  4.15.  Crosscorrelation  metrics  for  a  piecewise  basis  with  N  =  3, 
Gd  =  20  dB,  and  a  =  [12  4]^. 


Plkit) 

PCCL 

[dB] 

PoiW 

-17.394 

363.182 

P02it) 

-16.711 

363.182 

Pl2it) 

-14.061 

363.182 

and  4.16  have  been  calculated  from  these  correlation  functions  shown  in  Figures  4.20 
and  4.21  respectively. 
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Figure  4.21.  Crosscorrelation  functions  for  a  piecewise  basis  with 
AT  =  3,  Gd  =  40  dB,  and  a  =  [1  2  4]^. 


Table  4.16.  Crosscorrelation  metrics  for  a  piecewise  basis  with  iV  =  3, 
Gd  =  40  dB,  and  a  =  [1  2  Af. 


Plk{t) 

PCCL 

ldB| 

Poi{t) 

-26.603 

30061.120 

P02{t) 

-24.722 

30061.120 

Pl2it) 

-23.631 

30061.120 
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This  section  showed  that  the  piecewise  basis,  for  any  size  symbol  family,  is 
a  valid  basis.  This  basis  provides  slightly  better  performance  than  the  Chebyshev 
basis. 


4-4  Discrete  Code  Comparison 

Spread  spectrum  multiple  access  schemes,  such  as  Direct  Sequence  CDMA 
(DS-CDMA),  require  code  families  with  dispersive  crosscorrelation  properties  sim¬ 
ilar  to  those  used  in  NLS.  Unlike  the  Brown  symbols  described  earlier,  DS-CDMA 
utilizes  discrete  spreading  codes.  Gold  codes  were  invented  in  1967  specihcally  for 
use  in  multiple-access  spread  spectrum  systems  [30:135],  and  have  been  employed 
in  a  variety  of  applications  including  cellular  DS-CDMA.  Metrics  for  Gold  codes 
and  Simulated  Annealing  biphase  codes  are  provided  in  this  section  contrasting  the 
performance  of  Brown  symbols. 

The  codes  considered  in  this  section  are  hnite  in  time,  unlike  the  hnite  band¬ 
width  Brown  symbols.  The  metrics,  therefore,  have  been  calculated  for  discrete 
symbols  of  equivalent  signal  energy  and  approximately  equivalent  bandwidth  to  the 
Brown  symbols.  The  null-to-null  bandwidth  of  a  biphase  signal  is  approximately 
47r/rb,  where  Tb  is  the  bit  duration;  therefore  by  choosing  Tb  =  47r/Do  the  bandwidth 
is  approximately  that  of  the  Brown  symbols.  With  the  bit  duration  so  dehned,  the 
amplitude  of  the  biphase  symbol  is  then  chosen  so  that  the  energy  of  the  biphase 
symbol  is  equal  to  that  of  a  Brown  symbol. 

4-4D  Cold  and  SA  Biphase  Codes.  PRN-sequences  such  as  maximal-length 
codes  and  Gold  codes  are  commonly  employed  in  spread-spectrum  communications 
systems  [30:89].  Metrics  for  an  A^  =  3  family  of  31-bit  and  127-bit  Gold  codes  are 
illustrated  in  Tables  4.17  and  4.18  respectively. 
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Table  4.17.  Auto/Crosscorrelation  metrics  for  a  31-bit  Gold  code 
family  of  =  3. 


Plk{t) 

PSL 

ldB| 

ISL 

IdBI 

PCCL 

|dB] 

2 

^kk 

Poo(^) 

-17.786 

-7.227 

X 

13.857 

X 

pllif) 

-17.786 

-6.523 

X 

15.443 

X 

P22{t) 

-14.264 

-3.742 

X 

16.078 

X 

Poi{t) 

X 

X 

-9.827 

X 

44.600 

P02{t) 

X 

X 

-8.999 

X 

48.393 

Puit) 

X 

X 

-10.742 

X 

36.230 

Table  4.18.  Auto/Crosscorrelation  metrics  for  a  127-bit  Gold  code 
family  of  A^  =  3. 


Pikit) 

PSL 

[dB] 

ISL 

[dB] 

PCCL 

[dB] 

^kk 

^kl 

Poo(^) 

-24.014 

-6.023 

X 

267.233 

X 

Pii(f) 

-22.991 

-5.675 

X 

277.247 

X 

P22it) 

-19.154 

-2.240 

X 

325.511 

X 

Poi{t) 

X 

X 

-14.842 

X 

818.573 

P02{t) 

X 

X 

-15.228 

X 

701.845 

Pl2{t) 

X 

X 

-15.228 

X 

666.453 
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Codes  with  controlled  crosscorrelation  properties  have  recently  been  devised 
using  a  Simulated  Annealing  approach  [1:3-11].  Metrics  for  an  =  3  family  of 
31-bit  and  127-bit  SA  codes  are  illustrated  in  Tables  4.19  and  4.20  respectively. 

Table  4.19.  Auto/Crosscorrelation  metrics  for  a  31-bit  SA  code  fam¬ 
ily  of  A^  =  3. 


Plk{t) 

PSL 

[dB] 

ISL 

[dB] 

PCCL 

|dB] 

2 

^kk 

(^kl 

Poo{i) 

-14.264 

-3.106 

X 

19.472 

X 

Puii) 

-14.264 

-1.946 

X 

19.829 

X 

P22{t) 

-14.264 

-2.820 

X 

18.038 

X 

Polit) 

X 

X 

-11.765 

X 

41.693 

P02{t) 

X 

X 

-11.765 

X 

46.396 

Pl2{t) 

X 

X 

-11.765 

X 

51.002 

Table  4.20.  Auto/Crosscorrelation  metrics  for  a  127-bit  SA  code  fam¬ 
ily  of  A^  =  3. 


Pikit) 

PSL 

[dB] 

ISL 

[dB] 

PCCL 

[dB] 

2 

^kk 

(^kl 

Pmit) 

-19.154 

-2.793 

X 

356.686 

X 

PiiW 

-19.154 

-1.844 

X 

321.849 

X 

P22{t) 

-19.154 

-1.388 

X 

351.224 

X 

Poi{t) 

X 

X 

-16.501 

X 

763.032 

P02{t) 

X 

X 

-16.501 

X 

784.979 

Pl2{t) 

X 

X 

-16.501 

X 

778.029 

This  section  provided  comparative  results  for  biphase  codes  of  equivalent  signal 
energy  and  approximate  bandwidth.  Even  127-bit  Gold  and  SA  codes  do  not  perform 
as  well  as  Brown  symbols  in  terms  of  PCCL  and  PSL  indicating  that  Brown  symbols 
are  likely  to  be  well  suited  to  NLS  applications. 


4-5  Summary 

This  chapter  synthesized  and  analyzed  Brown  symbols,  using  the  process  out¬ 
lined  in  Section  3.7.2.  The  analytical  approach  described  in  Appendix  B  lead  to 
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the  construction  and  simulation  of  various  Brown  symbols  in  Matlab®.  Section  4.2 
considered  the  design  of  each  of  the  components  to  a  Brown  symbol:  the  Hermit 
Problem  solution,  phase-rate  basis  function  design,  phase  function  calculation,  and 
symbol  synthesis.  Section  4.3  provided  results  indicating  the  insensitivity  of  the  per¬ 
formance  of  Brown  symbols  to  the  choice  of  phase-rate  basis  functions  and  envelope, 
and  also  considered  physical  realizability.  These  results  support  the  claim  that  there 
is  no  sensitivity  in  rms  time  duration  of  the  correlation  functions  to  the  choice  of 
phase-rate  basis  functions,  so  long  as  the  criteria  of  optimality  and  orthonormality 
are  met.  Finally,  Section  4.4  contrasted  the  performance  of  Brown  symbols  with 
that  of  extant  biphase  symbols  with  Brown  symbols  showing  superior  performance 
in  rms  time  duration  characteristics  of  the  correlation  functions. 

Tables  4.21  and  4.22  provide  a  summary  of  the  comparative  performance.  As 
shown  in  Table  4.21,  the  autocorrelation  function  that  provides  the  best  overall 
performance  corresponds  to  the  Brown  symbol  family  with  the  polynomial  basis, 
which  was  derived  with  the  optimal  envelope.  The  lowest  PCCL  for  each  symbol 
family  is  extracted  and  displayed  in  Table  4.22.  The  Brown  symbols  derived  from 
the  non-ideal  uniform  taper,  and  Chebyshev  and  piecewise  basis  provide  the  best 
results.  Notably,  Brown  symbols  show  greater  mutual  dispersion  than  127-bit  Gold 
and  SA  codes  suggesting  that  they  are  better  suited  to  NFS. 
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Table  4.21.  Autocorrelation  metrics  for  various  symbol/code  fami¬ 
lies. 


Symbol  Basis/ 
Code 

PSL 

|dB] 

ISL 

|dB] 

^kk 

Sinusoidal 

-23.000 

-18.931 

0.016 

Polynomial 

-23.048 

-18.975 

0.016 

Chebyshev 

-13.267 

-7.076 

316.033 

Piecewise 

-13.262 

-7.109 

63.203 

31-bit  Gold 

-17.786 

-7.227 

13.857 

127-bit  Gold 

-24.014 

-6.023 

267.233 

31-bit  SA 

-14.264 

-3.106 

19.472 

127-bit  SA 

-19.154 

-2.793 

356.6866 

Table  4.22.  Crosscorrelation  metrics  for  various  symbol/code  families 
with  N  =  3. 


Symbol  Basis/ 
Code 

PCCL 

[dB] 

Sinusoidal 

-14.925 

30000.016 

Polynomial 

-15.978 

30000.016 

Chebyshev 

-22.655 

30315.616 

Piecewise 

-23.631 

30061.120 

31-bit  Gold 

-8.999 

48.393 

127-bit  Gold 

-14.842 

818.573 

31-bit  SA 

-11.765 

51.002 

127-bit  SA 

-16.501 

784.979 
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V.  Conclusions  and  Recommendations 


5.1  Restatement  of  Research  Goal 

As  stated  in  Section  1.2,  the  goals  of  this  research  are  to: 

1.  Explain  the  formulation  of  Brown’s  theorem  [5]  whilst  introducing  a  simplihed 
linear  algebraic  notation; 

2.  Describe  a  systematic  process  for  designing  Brown  symbols  through  an  analysis 
of  Brown’s  work  and  by  theoretical  extension; 

3.  Provide  a  theoretical  analysis  for  the  digitization  of  Brown  symbols; 

4.  Evaluate  performance  of  Brown  symbols,  comprising  dispersive  gain,  through 
simulation  and  numerical  approximation;  and 

5.  Compare  the  performance  of  various  Brown  symbols  with  discrete  codes  through 
simulation. 

5.2  Conclusions 

5.2.1  Brown  Symbol  Sensitivity.  Optimal  Brown  symbols  require  phase 
rate  functions  to  be  described  by  a  linear  combination  of  basis  functions  that  are 
orthonormal  within  the  weighted  inner  product  space  incorporating  A^O,  dehned 
in  Eq  (3.17).  Chapter  IV  considered  two  bases  that  meet  this  requirement,  the 
sinusoidal  basis  described  in  Section  4.3.1  and  the  polynomial  basis  in  Section  4.3.2. 
By  relaxing  the  optimality  requirement  for  envelope  A{uj),  two  other  bases  were 
considered,  the  Chebyshev  polynomials  in  Section  4.3.3  and  piecewise  continuous 
basis  in  Section  4.3.4. 

The  sinusoidal  and  polynomial  bases  provide  a  means  of  generating  N  op¬ 
timally  mutually  dispersive  symbols,  with  dispersive  gain  providing  control  over 
crosscorrelation  dispersion.  Both  of  these  bases  achieve  the  Brown  theoretical  au¬ 
tocorrelation  and  crosscorrelation  performance  in  rms  time  duration,  demonstrating 
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the  validity  of  Brown’s  theorem.  Notably,  for  unity  dispersive  gain,  the  form  of 
the  crosscorrelation  functions  vary  despite  having  equivalent  rms  time  durations. 
The  sinusoidal  basis  produces  crosscorrelation  functions  with  a  strong  mainlobe  and 
dominant  sidelobes  at  approximately  -18  dB,  whereas  the  polynomial  basis  produces 
crosscorrelation  functions  with  a  broad  mainlobe  and  no  sidelobes.  When  the  disper¬ 
sive  gain  factor  is  increased,  however,  both  bases  produce  crosscorrelation  functions 
with  extensive  dispersion  at  approximately  -20  dB.  This  result  suggests  that  Brown 
symbols  need  only  be  considered  for  symbols  with  large  dispersive  gain,  which  also 
relates  to  large  time-bandwidth  products. 

Brown  symbols  show  insignihcant  degradation  due  to  windowing,  provided  the 
window  width  is  selected  according  to  dispersive  gain  as  illustrated  in  Figure  4.7. 
This  suggests  that  Brown  symbols  are  suitable  to  radar  applications.  When  com¬ 
pared  to  conjugate  LFM,  Brown  symbols  have  superior  autocorrelation  statistics  but 
produce  crosscorrelation  functions  with  a  consistent  10  dB  lower  PCCL.  Although 
the  LFM  pulse  pair  has  better  PCCL  performance.  Brown  codes  exhibit  superior 
performance  in  terms  of  rms  time  duration.  For  example,  the  comparison  of  pulses 
with  equivalent  rms  time  duration,  in  Figure  4.9,  showed  that  Brown  codes  exhibit 
greater  compression.  PCCL  is  an  excellent  indicator  of  a  symbol  families  suitability 
for  NLS,  however,  any  two  symbols  will  have  only  one  crosscorrelation  function  and 
can  claim  justly  to  be  mutually  dispersive;  this  is  not  so  for  larger  code  families. 
LFM  pulses  have  been  shown  to  provide  excellent  performance  in  NLS  but  are  of 
minimal  beneht  to  NLS  systems  due  to  the  severe  limitation  on  symbol  family  size, 
N  =  2.  Conversely,  Brown  symbols  are  scalable  to  any  sized  family,  they  perform 
better  in  terms  of  rms  time  duration  and  are  comparable  to  LFM  in  terms  of  PCCL. 

When  the  constraint  of  envelope  optimality  is  removed,  and  a  uniform  enve¬ 
lope  applied,  the  properties  of  the  autocorrelation  function  worsen.  This  is  shown 
in  Tables  4.9  and  4.13  with  a  signihcant  increase  in  rms  time  duration  and  a  10  dB 
increase  in  PSL.  The  Chebyshev  and  piecewise  bases,  due  to  the  increase  in  auto- 
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correlation  dispersion,  provide  better  crosscorrelation  dispersion  at  unity  gain  than 
both  the  sinusoidal  and  polynomial  bases.  The  dispersive  gain  can  be  increased  to 
improve  crosscorrelation  dispersion,  therefore,  small  differences  in  PCCL  are  easily 
accounted  for  in  this  way,  negating  the  requirement  to  use  a  non-ideal  envelope.  The 
shape  of  the  correlation  plots  is,  however,  signihcantly  different.  For  lower  dispersive 
gain,  or  time-bandwidth,  the  crosscorrelation  functions  resulting  from  these  bases  ex¬ 
hibit  much  flatter  dispersive  characteristics.  A  relaxing  of  the  optimality  constraint 
therefore  effects  both  the  autocorrelation  and  crosscorrelation  statistics  providing 
marginal  improvements  in  PCCL  and  dispersion,  potentially  flatter  crosscorrelation 
functions,  but  higher  sidelobes. 

5.2.2  Discrete  Code  Comparison.  Section  4.4  provided  a  brief  insight  into 
the  performance  of  extant  waveforms  that  are  suited  to  NLS.  Gold  codes  and  SA 
codes  were  considered  in  comparison  to  Brown  symbols.  Neither  the  31-  or  127-bit 
Gold  or  SA  codes  provide  comparable  performance  to  Brown  symbols.  Although  the 
127-bit  Gold  codes  have  comparable  PSL  values,  both  types  of  codes  exhibit  much 
worse  compression  and  dispersion  characteristics  than  Brown  symbols.  Also,  these 
codes  offer  no  means  of  controlling  the  degree  of  dispersion. 

5.3  Significant  Results  of  Research 

NLS  symbol  families  constructed  from  signals  with  desirable  correlation  prop¬ 
erties,  such  as  LFM  symbols  [28]  and  discrete  Gold  codes  [2],  have  already  been 
implemented.  However,  LFM  pulse  families  are  severely  limited  in  size  and  the  per¬ 
formance  of  Gold  codes  is  constrained.  Theoretically,  Brown  symbols  can  be  easily 
tailored  for  utilization  in  NLS.  The  theorem  claims  that  these  symbols  are  optimally 
mutually  dispersive,  but  has  until  now  remained  untested. 

This  research  reintroduced  Brown’s  theorem  relating  the  system  model  as¬ 
sumptions  and  theory  formulation.  As  a  result  of  this  research,  a  systematic  method 
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for  synthesizing  Brown  symbols  is  now  defined  and  demonstrated.  The  synthesis 
and  simulation  process  provides  the  theoretical  background  necessary  to  enable  dig¬ 
itization  of  Brown  symbols.  The  Brown  symbols  that  were  created  and  analyzed 
achieve  theoretical  performance  and,  for  large  dispersive  gain,  exhibit  desirable  cor¬ 
relation  characteristics.  The  dispersive  gain  factor  was  introduced  and  integrated 
into  the  formulation  and  design  of  Brown  symbols,  thus  providing  a  crosscorrelation 
dispersion  control  parameter.  Subsequently,  Brown  symbols  have  been  shown  to  be 
well-suited  for  physical  implementation  into  pulse  Doppler  radar  systems,  exhibiting 
more  desirable  correlation  statistics  than  comparable  discrete  codes.  As  such,  the 
goals  stated  in  Section  1.2  have  been  accomplished. 

5.4  Recommendations  for  Future  Research 

Recommendations  for  future  work  center  on  expanding  the  complexity  of  the 
radar  model  and  further  demonstrating  Brown  symbol  performance,  robustness  and 
sensitivity.  The  Brown  theorem  can  also  be  expanded  to  image  processing. 

The  radar  model  referenced  in  this  research  incorporated  simplifying  assump¬ 
tions  of  a  noiseless,  distortion  free  environment  with  stationary  single-point  targets. 
To  better  understand  the  suitability  of  Brown  symbols  to  pulse  Doppler  radar  these 
assumptions  should  be  altered  to  better  reflect  a  real  radar  environment.  As  such, 
future  work  should  consider  the  implementation  of  Brown  symbols,  with  different 
generating  bases,  in  the  presence  of  noise.  Also,  Brown  symbols  appear  reliant  on 
the  integrity  of  the  phase  functions,  therefore,  an  analysis  of  the  effect  of  Doppler 
and  clutter  on  the  properties  of  the  correlation  functions  would  be  desirable. 

Brown  [5:13]  has  proposed  a  study  of  Brown  symbols  for  implementation  with 
functions  of  more  than  one  variable,  such  as  images.  For  example,  future  work  could 
consider  expanding  the  theory  of  mutually  dispersive  symbols  to  2-dimensions  and 
analyze  the  performance  of  signals  generated  in  that  way.  A  possible  application 
identified  by  Brown  is  the  coded-multiplexing  of  color  television  signals  onto  a  single 
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channel.  Other  2-diniensional  research  could  include  the  expansion  of  the  theory  to 
facilitate  the  design  of  radially  symmetric  Brown  surfaces. 

The  principle  of  stationary  phase  may  be  employed  to  approximate  the  inverse 
Fourier  transform  of  the  crosscorrelation  functions.  By  considering  this  approach,  an 
intuitive  avenue  of  research  becomes  apparent,  i.e.,  the  approximation  to  the  cross¬ 
correlation  function  is  a  linear  combination  of  complex  sinusoids;  their  phases  are 
determined  by  the  stationary  points  of  the  phase-rate  function.  Periodic  phase-rate 
functions,  such  as  the  sinusoidal  basis,  are  therefore  likely  to  generate  crosscorre¬ 
lation  functions  with  constructive  superposition  leading  to  sidelobes.  A  random 
phase-rate  process  with  uniform  distribution  would  considerably  alleviate  this  issue. 
An  expansion  to  the  theory  provided  in  this  thesis  should  consider  constructing  a 
framework  for  the  use  of  random  processes  as  basis  functions.  Initial  work  in  this 
area  shows  crosscorrelation  functions  with  excellent  dispersive  properties  and  PCCLs 
on  the  order  of  -32  dB. 
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Appendix  A.  Ambiguity  Function  Analysis 

A.l  Ambiguity  Surfaces 

This  appendix  displays  a  series  of  ambiguity  surfaces  generated  for  various 
radar  waveform  types. 

A.  1.1  Single  Pulse  Waveform.  The  most  basic  pulsed  radar  consists  of  a 
single  pulse  modulated  onto  a  carrier. 

x(t)  =  rect  (A.l) 

The  ambiguity  surface  for  this  function,  with  /c  =  0,  is  determined  from 

X(r,  A  =  J  rect  rect  (A.2) 

and  is  illustrated  in  £g.  A.l. 


A. 1.2  Pulse-Train  Waveform.  The  next  most  advanced  pulsed  radar  con¬ 
sists  of  a  train  of  pulses  modulated  onto  a  carrier,  used  to  improved  detection  and 
Doppler  resolution.  In  practise  an  infinite  pulse-train  is  not  feasible,  therefore  in  this 
analysis  an  N  pulse  train  is  considered. 

x{t)  =  ^  rect  (- — Qp'^fot  (A. 3) 


The  ambiguity  surface  for  this  function,  with  /c  =  0,  is  determined  from 


x(d  = 


rtt  -<=‘  (  ( A-4) 


and  is  illustrated  in  fig.  A.2  for  N  =  6. 
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Figure  A.l.  Ambiguity  surface  of  a  single  pulse  radar  waveform. 


Figure  A. 2.  Ambiguity  surface  of  a  pulsed  radar  waveform. 
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A.  1.3  Linear  FM  Waveform.  The  linear  FM  pulse  also  has  interesting 
properties,  as  discussed  in  section  2. 3. 3. 2.  This  signal  sweeps  a  frequency  range,  <F, 
and  can  be  described  by 


a^LFM(^)  =  rect 


where  fc  =  $/T 


(A.5) 


The  ambiguity  surface  for  this  function  is  determined  from 

x(r, ")  =  y  ‘■ect  A  jj" ““rect 


(A.6) 


and  is  illustrated  in  £g.  A. 3. 


Figure  A. 3.  Ambiguity  surface  of  a  single  pulse  LFM  waveform. 
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Appendix  B.  Basis  Function  Derivations 

B.l  Sinusoidal  Basis 

The  sinusoidal  functions  are  an  obvious  choice  to  form  a  bases  due  to  their 
orthogonality  properties.  With  the  optimal  envelope  for  Brown  symbols,  however, 
the  orthogonality  must  hold  within  the  weighted  inner  product  space,  (dfc,  'dz)^2.  The 
weighted  inner  product  of  sinusoidal  functions  is  dehned  such  that 

dk{oj)  =  6sin(afcC4;),  (B.l) 


then 


{dk,di)A2  =  7^  I  dk{uj)di{u)A^{uj)duj 

ZTT 


1 

Qq 


^0 

[  b  sm{ak^)b  sin(azci;)  —  cos^ 
-^0  i  Zn 

2 

62 

/  sin(aA:Cc;)  sin(a«a;)(ia;  +  — 

'  -^0  iZ, 

2 


£)*• 


2 


0 


While  afc  are  integers  for  all  k,  and  a^  7^  a^.  Thus,  the  sinusoidal  functions  in  the  hrst 
term  of  Eq  (B.2)  are  orthogonal  when  evaluated  over  an  integer  number  of  cycles. 
Subsequently,  a  sufficient  constraint  for  orthogonality  here  is  that  hlo  =  2m7r,  for 
any  integer  m.  The  weighted  norm,  k  =  I,  reduces  the  first  term  as  such 


Slid  {akuj)du 


F 

2ffo 


sin(afcCi;) 


1  2 


2  ■ 


(B.3) 


The  second  term  is  simply  an  integral  of  a  cosine  function  over  one  complete  cycle, 
which  will  always  be  zero.  Therefore  the  sinusoidal  basis  functions  are  orthogonal. 
By  choosing  b  so  that  Eq  (B.3)  equals  unity,  an  orthonormal  basis  is  described. 
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Therefore  an  orthonormal  sinusoidal  basis  is  dehned  by 


where 


dfc(a;)  =  A/2sin(afca;), 


hlo  =  217171 

a.k,  m  G  Z 

7^  a;,  V  /c  ^ 


(B.4) 


(B.5) 


B.2  Polynomial  Basis 

There  are  a  number  of  orthogonal  polynomial  bases  such  as  Chebyshev  func¬ 
tions,  Legendre  functions,  Hermite  functions  and  so  on.  The  suitability  of  each  of 
these  alternate  bases  is  dependent  on  them  being  odd  functions  that  retain  their 
orthogonality  within  the  weighted  inner  product  space,  The  search  for 

a  suitable  basis  begins  by  dehning  the  basis  element  as  a  weighted  set  of  odd 
polynomial  functions  of  order  2p  —  1. 

p-i 

M^)  =  (B.6) 

n=0 

then 


271 


d  fc  (o; ) -dz  (o; )  (ci; )  dcu 


1 

Ytt  /- 


—  OO 

p— 1 


J  ^  ^  ^  ^  cos^  ( '^  ]  doj 


2 

hlo 


2  n=0 

p—1  p— 1 


p  ^  p  ^  p 

EE  bfcnb/m  / 

n=0  m=0  ~ 


m=0 

yo 
2 


a 


\n 


cos' f  ^  )  dcu.  (B.7) 

-^0  \  “0 
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The  resultant  of  the  integral  is  a  function  only  of  m  and  n  and  is  therefore  redehned 


^0  -  , 

Emn  =  -^  ,;2(n+„^+l)  ^^^2  (^]du;,  (B.8) 

VloJz^  V^o/ 

which  simplihes  Eq  (B.7)  to 

p—1  p—1 

EE  (B,9) 

n=0  m=0 

Again  this  may  be  further  simplihed  using  linear  algebraic  notation  as 

{^k,^i)A^  =  hlEhi.  (B.IO) 

Dehning  B  =  [bg  bi  b2  . . .  b^r-i]  permits  further  simplihcation  of  Eq  (B.IO).  Solving 
for  ort honor mality  provides 

=  B'^EB 

=  liv.  (B.ll) 

Therefore  when  E  is  invertible 


BB^  =  E-\  (B.12) 

or  in  other  words  with  B^  as  the  Cholesky  decomposition  matrix  to  the  inverse  of 
E  another  orthonormal  basis  within  the  weighted  inner  product  space  is  dehned. 

B.3  Chebyshev  Functions 

Changing  the  spectral  envelope,  A{u),  introduces  a  different  weighted  inner 
product  space  within  which  to  design  an  orthonormal  basis.  An  envelope  that  is 
uniform  across  the  spectral  bandwidth  invites  the  use  of  Chebyshev  polynomials. 
This  is  because  any  pair  of  Chebyshev  polynomials  exhibit  the  orthogonality  prop- 
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erties  described  in  Eq  (B.13)  [36:158]: 


Tm{uj)Tn{uj) 


duj 


0,  m  ^  n 

■*  TT,  m  =  n  =  0 

7r/2,  otherwise 


(B.13) 


where  Tk{uj)  is  an  order  Chebyshev  polynomial.  The  odd  Chebyshev  functions 
can  be  dehned  as 


T2m\-l{oj)  =  COs[(2m  +  1)  COS 

fc=0  ' 

With  A{uj)  dehned  as  a  uniform  envelope  on  the  interval  [—r2o/2,  flo/2]  as  in  Eq 
(B.15),  the  Chebyshev  functions  are  scaled  to  induce  orthonormality  within  this 
inner  product  space  as  shown  in  Eq  (B.16). 


A{uj)  = 

ldn{uj)  = 


(B.15) 

(B.16) 


Therefore  the  weighted  inner  product  of  any  pair  of  basis  functions  is 


(■dfc,  'd«)A2 


1 

1 


dfe  (a;)'dz  {uj)A‘^  {uj)du 


— du 


(B.17) 
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With  a  change  of  variable  (7  =  2uj/VLq)  the  integral  takes  on  the  familiar  form  of 
Eq  (B.13),  where  for  all  k,  I, 


2  f  T2k+i{'y)T2i+i{'y) 


Ski- 


(B.18) 


Therefore  the  Chebyshev  basis  as  dehned  in  Eq  (B.16)  is  orthonormal  for  a  non-ideal 
uniform  taper. 


B.4  Piecewise  Continuous  Basis 

LEM  pulses  were  originally  demonstrated  in  the  early  implementation  of  a  NLS 
scheme.  These  pulses  have  linear  phase-rate  functions  with  spectral  content  that 
is  approximately  uniformly  dispersed  over  a  compact  frequency  bandwidth.  The 
duality  of  the  Fourier  transform  provides  that  functions  of  frequency  having  linear 
phase-rate  functions  will  correspond  to  time  signals  that  are  uniformly  dispersed  and 
compactly  supported.  This  is  a  desirable  property  of  the  crosscorrelation  functions 
as  the  presence  of  sidelobes  is  expected  to  be  signihcantly  mitigated.  This  section 
explores  the  use  of  piecewise  linear  phase-rate  functions,  such  as  the  saw-tooth  wave¬ 
form.  To  provide  a  means  of  analysis  this  basis  utilizes  the  same  uniform  taper  as 
the  Chebyshev  basis  as  described  in  B.15.  As  such,  the  function  that  describes  the 
normalized  basis  is  dehned  as 


'&k{u!)  = 


4 


f  I  n  /o^  ^  +  ^o/‘^  —  ^k/8 

(a,  +  n„/2)rect(^ - - 


2aj.  — 1 

^  (-l)”(a;  -f  a;n)rect' 


n=l 


UJ  —  Ur, 

U/2 


+  (a,  -  SJ„/2)rect(^ - - 


(B.19) 


where 


ktk  —  f^o/ Hk 
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y(n/afc  -  1). 


That  is,  afc  is  the  number  of  cycles  over  hlo;  is  the  period,  and  Un  are  the  zero 
crossings  as  illustrated  in  Figure  B.4. 


Figure  B.l.  Saw-tooth  basis  function  with  a*,  =  4  and 

fio  =  2. 

The  selection  of  is  constrained  to  ensure  that  ^k{^)  forms  an  orthonormal 
basis.  This  can  be  accomplished  by  ensuring  that  max  {a?;,  aj}/min{afc,  a^}  is  radix 
2  for  all  k,  1.  Figure  B.4  illustrates  that  the  product  of  two  basis  functions  with 
max  {afc,  a;}/min{afc,  a;}  =  2  produce  equivalent  positive  and  negative  areas. 
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Figure  B.2.  Orthogonal  piecewise  functions. 

'(9fe(cc;)  and  selected  such  that 

maxjafc,  az}/min{afc,  a;}  =  2,  (a),  and 
Area  under  (b). 
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